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ABSTRACT 


Ocean  acoustic  tomography  permits  the  mapping  of  various  properties  of  a  body 
of  water  through  indirect  means.  The  technique  utilizes  travel  time  variations  for  an 
acoustic  signal  to  determine  the  structure  of  the  ocean  medium  via  inverse  mathematical 
methods.  The  scale  of  any  tomography  experiment  is  fundamentally  limited  by  the 
signal  to  noise  ratio  at  the  receiver.  Through  the  use  of  a  near  vertical  acoustic  array, 
normal  mode  modeling  of  the  Leal  environment  and  a  modal  beamformer,  array  gains 
are  possible  which  greatly  extend  the  maximum  separation  between  source  and  receiver. 
Additionally,  the  technique  provides  temporal  resolution  of  the  modal  components  of 
the  arriving  signal. 

A  time  domain  modal  beamformer  for  a  near  vertical  acoustic  arvay  has  been 
developed.  It  has  realized  a  nominal  a'ray  gain  of  6  dB  for  the  Heard  Island 
Experiment  vertical  array  deployed  off  California.  The  primary  obstacle  to  the 
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I .  INTRODUCTION 


A.  THESIS  SUMMARY 

The  objective  of  this  thesis  is  to  develop  a  software 
package  to  beamform  acoustic  signals  used  in  ocean  acoustic 
tomography.  The  goal  of  tomography  signal  processing  is  the 
precise  measurement  of  acoustic  travel  time.  Sound  speed  is 
a  well  understood  function  of  temperature,  pressure  and 
salinity.  Therefore,  a  fluctuation  in  acoustic  travel  time  is 
indicative  of  changes  in  the  environment  through  which  the 
acoustic  energy  has  passed.  Once  the  travel  time  fluctuations 
are  known,  inverse  mathematical  techniques  can  be  used  to 
infer  various  properties  of  the  ocean  medium.  [Ref.  1] 

One  source  of  travel  time  measurements  is  through  the  use 
of  explosive  devices.  Although  such  signals  are  easy  to 
generate,  they  may  not  provide  the  required  precision  owing  to 
dispersion  of  the  various  frequency  components  in  the 
impulsive  signal.  Additionally,  such  signals  are  difficult  or 
impossible  to  replicate.  [Ref.  2] 

A  better  method  that  has  been  employed  in  recent  years  is 
the  use  of  maximal -length  sequences.  The  technique  utilizes 
a  pseudorandom,  binary  sequence  as  the  phase  modulation  for  an 
electronically  generated  bandpass  signal.  Maximal-length 
sequences  are  well  suited  to  travel  time  measurements  by 
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virtue  of  their  deterministic  nature,  autocorrelation 
properties  and  simplicity. 

The  goal  of  the  work  described  in  the  following  pages  is 
to  provide  sufficient  gain  to  permit  the  recovery  and  decoding 
necessary  for  accurate  travel  time  measurements  for  acoustic 
paths  of  extreme  length  (in  this  case  10,000  nautical  miles). 
Specifically,  the  programming  package  should  be  able  to: 


•  Exploit  the  spatial  structure  of  the  incident  wave  field 
to  discriminate  among  signals  arriving  from  different 
directions . 

•  Utilize  information  regarding  time  varying  array  tilt  and 
depth  to  estimate  the  array  geometry. 

•  Provide  a  stable  virtual  array  by  using  the  array  geometry 
and  modal  structure  of  the  immediate  environment. 

•  Distinguish  among  the  various  modal  components  of  the 
target  signal. 

•  Provide  both  time  domain  and  frequency  domain  analysis. 

•  Remain  sufficiently  flexible  so  as  to  permit  the 
processing  of  arrays  of  differing  construction  in  future 
tomography  experiments. 

•  Provide  some  measure  of  fault  tolerance  with  respect  to 
the  receiving  array. 


The  remainder  of  this  thesis  is  structured  as  follows: 


•  Chapter  II. 

•  Chapter  III. 

•  Chapter  IV. 

•  Chapter  V. 


Acoustic  Wave  Propagation  Theory 
Modal  Beamforming 
The  Heard  Island  Array 
Results  and  Conclusions 
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Chapter  II  presents  an  introduction  to  acoustic  wave 
propagation.  The  approach  taken  is  that  of  normal  mode 
propagation  in  a  range  independent  channel. 

The  third  chapter  is  an  introduction  to  the  concept  of 
modal  beamforming.  It  reviews  frequency  domain  and  time 
domain  beamformers.  The  algorithm  utilized  in  this  study  is 
derived . 

Chapter  IV  describes  the  construction  and  subsequent 
deployment  of  the  receiving  array  used  in  this  experiment. 
Data  acquisition  and  preprocessing  of  the  acoustic  signal  are 
discussed . 

The  final  chapter  presents  a  description  of  the  array 
dynamics  encountered  and  frequency  domain  output  from  the 
software  developed.  Additionally,  proposals  for  system 
integration  and  performance  improvements  are  detailed. 

Appendix  A  contains  the  source  code  which  has  served  as 
the  body  of  this  work.  In  addition  to  the  beamformer,  various 
utility  programs  are  included  which  provide  for  array  geometry 
data  reduction. 

B.  THE  HEARD  ISLAND  EXPERIMENT 

The  performance  of  the  software  package  is  evaluated  on  a 
raw  data  set  (vice  synthetic  data).  The  data  was  acquired 
during  the  Heard  Island  Experiment  which  occurred  during  the 
period  January  23  -  February  2,  1991.  The  purpose  of  the 
experiment  was  to  determine  the  reliability  of  global  acoustic 
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paths  for  tomographic  analysis.  Specifically,  it  is  desired 
to  establish  the  viability  of  these  transmission  paths  for  a 
proposed  multi-national  attempt  to  detect  a  decreasing  trend 
in  acoustic  travel  time.  Such  a  change  would  indicate  an 
overall  warming  along  the  path,  providing  the  first  convincing 
evidence  to  the  existence  of  global  warming  [Ref.  1]. 

The  signals  were  transmitted  from  the  vicinity  of  Heard 
Island  in  the  southern  Indian  Ocean.  The  site  was  selected 
because  it  is  central  to  a  web  of  open  water  paths  extending 
through  all  the  worlds  oceans  as  shown  in  Figure  1.1. 


The  transmitting  ship  was  the  R/V  Cory  Chouest.  The 
transmitter  utilized  a  ten  element  vertical  array  with  a 
maximum  of  five  elements  active  at  any  time.  A  nominal  source 
level  of  209  dB  re  1  pPa  at  1  meter  was  realized.  Individual 
projectors  employed  were  HLF4LL  very  low  frequency  sources. 
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A  variety  of  signal  types  were  sent,  including  continuous 
wave  and  maximal -length  sequences.  A  carrier  frequency  of  57 
Hz  was  chosen  both  for  its  low  absorption  losses  and  the 
ability  to  distinguish  it  from  the  50  and  60  Hz  frequencies 
generated  by  power  plants  world  wide. 

The  Monterey  component  of  the  Heard  Island  Experiment 
involved  a  collaboration  between  the  Naval  Postgraduate 
School,  the  Monterey  Bay  Aquarium  Research  Institute,  the 
Massachusetts  Institute  of  Technology,  Woods  Hole 
Oceanographic  Institution,  SAIC  and  the  Moss  Landing  Marine 
Laboratory.  The  R/V  Point  Sur  deployed  a  32  element  vertical 
array  approximately  70  nautical  miles  south-west  of  Monterey 
Bay . 
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II.  ACOUSTIC  WAVE  PROPAGATION  THEORY 


A.  ACOUSTIC  PROPAGATION  IN  AN  OCEAN  WAVEGUIDE 
1 .  The  Inhomogeneous  Wave  Equation 

The  homogeneous  linear  wave  equation  [Ref.  3], 


V»p  =  —  , 

dt^ 


(2.1) 


incorporates  two  implicit  assumptions.  The  first  restricts 
the  sound  speed  to  a  constant  value  in  the  region  of  interest. 
Following  the  development  of  Coppens  [Ref.  4],  if  the  sound 
speed  is  now  permitted  to  vary  in  space  (but  not  time),  then 
c=c(x,y,z)  may  be  substituted  into  the  wave  equation  without 
loss  of  generality. 

The  second  implicit  assumption  restricts  the  use  of 
the  wave  equation  to  sound  fields  which  are  free  of  sources. 
Before  including  a  source  term,  the  selection  of  a  specific 
coordinate  system  is  appropriate. 

Given  that  sound  in  the  ocean  does  not  spread 
spherically  in  a  free  field,  but  radially  with  upper  and  lower 
boundaries,  implies  the  use  of  a  cylindrical  coordinate 
system.  To  further  simplify  the  problem  to  one  sufficient  for 
this  work,  the  sound  field  shall  be  assumed  to  exhibit  radial 
symmetry  about  the  source,  and  sound  speed  variations  shall  be 
restricted  to  depth,  c(z).  Based  on  these  restrictions. 
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application  of  the  Lapacian  yields  the  homogeneous,  range 
independent  wave  equation, 


1  d  I  d  ^  _  1  ^ 

r  dr[  drr  c^(z) 


8^ 


(2.2) 


(z)  dt^' 

The  inclusion  of  a  source  term  requires  that  the 
inhomogeneous  wave  equation  reduce  to  homogeneous  form  in 
regions  which  are  free  of  sources  of  acoustic  energy.  This 
requirement  implies  that  the  delta  function  for  a  point  source 
located  at  r=0  and  z=z^  have  the  form  [Ref.  3] 


6(r-r^)  =  ^^5(r)  6(z-z^)  ,  (2.3) 

2 11  r 

where  r  is  the  spherical  radius  vector,  r  is  the  cylindrical 
radial  coordinate  and  z  is  the  cylindrical  depth  coordinate. 

Inclusion  of  this  term  in  the  wave  equation  yields  the 
Helmholtz  equation  for  a  monofrequency  point  source  of  unity 
amplitude , 


’i  a 

r  — 1 

.  J  ^  f 

r  dr 

dr] 

az"  \c(z)j 

--biz)  biz-zj 
r 


(2.4) 

2.  Solutions  to  the  Inhomogeneous  Wave  Equation 

Having  adopted  a  cylindrical  coordinate  system, 
equation  (2.4)  can  be  solved  by  separation  of  variables. 
Furthermore,  the  complete  solution  can  be  treated  as  a  linear 
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combination  of  normal  modes, 


p{r.  z.  t)  =  Rjz)  ZJz)  .  (2.5) 

in 

The  normal  modes  (Z^^)  form  an  orthonormal  set  of 
eigenfunctions  in  2  which  satisfy  the  sourceless  Helmholtz 
equation, 


dz^ 


(2.6) 


subject  to  the  appropriate  surface  and  bottom  boundary 
conditions,  and  the  orthonormal  condition, 

fzjz)  Z^,(z)  dz  =  6^,  (2.7) 

where  is  the  eigenvalue  for  the  m^*^  eigenfunction  (or  normal 
mode)  and  6„  is  the  Kronecker  delta  function . 

A  closed  form  expression  for  the  R„(r)  may  be  obtained 
by  substitution  of  equation  (2.6)  into  the  inhomogeneous  wave 
equation  (2.4),  multiplication  of  all  terms  by  Z^,  integrating 
over  2  and  application  of  the  normalization  condition  (2.7). 
The  result,  after  manipulation,  is 


1  _d 
I  dr 


-l6(r)Z„(z,), 


(2.8) 


the  inhomogeneous  Bessel's  equation.  It  has  the  known 
solution 

R^ir)  =  -7  7iZ„(2^) (K^r)  ,  (2.9) 
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where  is  the  Hankel  function  of  the  second  kind  and 
order  zero.  Substitution  of  this  form  into  equation  (2.5) 
yields 

p{r.z,t)  =  (K^r)  .  (2.10) 

m 

Closed  form  solutions  for  the  Zj,(z)  are  either 
difficult  or  impossible  for  all  but  select  cases. 
Fortunately,  efficient  numerical  algorithms  exist  which 
provide  for  the  rapid  solution  of  the  depth  dependent 
functions.  [Ref.  5] 
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III.  MODAL  BEAMFORMING 


A.  PRELIMINARY  CONCEPTS 

Propagating  waves  can  be  modeled  as  functions  of  both 
space  and  time.  Consequently,  the  attributes  of  such  models 
can  be  used  to  extract  information  from  real  wave  fields. 
Beamforming  exploits  the  temporal  and  spatial  characteristics 
of  a  specified  environment  to  enhance  a  particular  aspect  of 
the  wave  field  while  attenuating  undesirable  components.  Such 
an  operation  can  be  loosely  termed  constructive  reenforcement 
of  the  desired  signal  or  noise  rejection.  The  information 
exploited  in  this  development  includes  the  modal  structure  of 
the  immediate  environment  and  the  instrument  data  supplied  by 
the  array  itself. 

B.  TIME  DOMAIN  BEAMFORMING 

The  acoustic  signal  received  at  a  given  hydrophone  can  be 
expressed  as 

p{t,x,y,z)  =  S{t,x,y,  z)  *N{t,x,y,  z)  ,  (3.1) 

where  p  is  the  pressure  at  the  hydrophone,  S  is  the  desired 
signal  and  N  is  the  local  noise  field.  All  are  functions  of 
time  and  hydrophone  location. 

If  one  now  considers  the  signal  received  at  an  array  of 
hydrophones,  equation  (3.1)  takes  the  form 


10 


P(t)  =  Sit)  *Nit)  , 


(3.2) 


where  p,  S  and  N  are  all  N  x  l  vectors  representative  of  the 
individual  array  elements.  The  location  of  individual 
elements  is  implied  by  the  vector  index  for  a  given 
hydrophone . 

Classical  (or  plane  wave)  beamforming  incorporates  the  use 
of  a  complex  steering  vector  (E)  to  impose  a  phase  shift  on 
the  individual  elements  in  order  to  enhance  the  sensitivity  of 
the  array  to  signals  propagating  in  a  specific  direction, 


E  = 


(3.3) 


The  output  signal  from  such  an  array  is 

bfp^it)  =  E^iSit)  -^Nit))  .  (3.4) 

The  subscript  (pw)  indicates  that  the  beamformer  assumes  the 
presence  of  plane  waves.  Additionally,  the  implicit 
assumption  is  made  in  the  above  expression  that  the 
autocorrelation  function  of  the  noise  field  will  be 
identically  zero  for  any  two  (distinct)  hydrophones.  The 
success  or  failure  of  the  beamformer  will  be  a  direct 
reflection  of  the  validity  of  these  assumptions. 
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If  the  desired  signal  is  not  monofrequency,  but  of  finite 
bandwidth,  simple  phase  shifting  may  not  be  suitable.  This  is 
particularly  true  if  the  target  signal  is  carrying  information 
such  that  the  communication  bit  time  is  comparable  to,  or 
greater  than  the  time  taken  for  the  wave  to  propagate  across 
the  array.  Under  such  circumstances,  the  beamformer  would 
impose  phase  shifts  in  distinct  communication  bit  segments 
simultaneously.  The  resulting  distortion  to  the  bandpass 
signal  is  unacceptable  in  most  communications  applications. 
Beamforming  under  these  conditions  requires  the  application  of 
true  time  delays  vice  phase  shifts. 

Assuming  an  array  of  spatial  extent  such  that  time  domain 
beamforming  is  required  implies  a  summation  of  the  form 

N 

=^{S(t-Xj.)*N{t-x^)).  (3.5) 

where  n  is  the  index  used  to  describe  the  array  elements  and 
represents  the  steering  delay  applied  to  the  n^^  element 
[Ref.  6].  This  arrangement  places  no  additional  restrictions 
on  the  target  signal  characteristics  or  array  construction  in 
a  given  physical  system. 

C.  THE  MODAL  BEAMFORMER 

The  signal  model  employed  above  is  incomplete,  in  that  it 
does  not  account  for  the  nonuniform  distribution  of  acoustic 
energy  in  the  sound  channel.  Vertical  arrays  inherently 
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sample  the  modal  structure  of  the  immediate  environment. 
Therefore,  one  needs  to  take  the  information  provided  by 
normal  mode  theory  into  account  when  the  acoustic  array  spans 
a  region  of  nonuniform  energy  density. 

The  appropriate  refinement  to  the  complex  steering  vector 
which  incorporates  the  modal  nature  of  the  signal  arrival  is 


.  .  Z„(z2)e''*“ 

where  the  m^^  column  of  the  matrix  is  the  generalized  steering 
vector  for  the  m^^  normal  mode.  The  individual  steering 
vectors  now  include  amplitude  factors  (Z^(z^))  that  reflect  the 
value  of  the  eigenfunction  at  the  various  hydrophone  depths 
[Ref.  7].  The  n^^  row  contains  the  values  of  those  steering 
vectors  at  the  depth  of  the  n^*’  hydrophone.  The  phase  terms 
account  for  phase  shifts  required  for  the  beam  steering 
applicable  to  all  modes  (as  required  by  array  geometry)  and 
phase  delays  among  the  various  modes  due  to  their  individual 
propagation  speeds.  Specifically,  the  take  the  form 

Ur„„  =  0„+a„,  (3 .1) 

where  6^  is  determined  by  the  hydrophone  location  and  is 
determined  by  the  differences  in  propagation  speeds  over  the 
transmission  path  for  the  normal  modes  present. 


Modifying  the  steering  matrix  to  implement  time  delays 
(vice  phase  shifts)  is  straight  forward.  A  beamformer  could 
be  implemented  by  summing  the  output  of  individual  hydrophones 
(p(t))  at  this  point.  However,  to  do  so  implies  a 
deterministic  evaluation  of  the  a^.  Such  an  evaluation  would 
require  near  perfect  knowledge  of  the  conditions  over  the 
transmission  path  prior  to  commencement  of  the  experiment. 
Given  that  the  aim  of  tomography  is  the  mapping  of  these 
conditions,  such  an  approach  is  inappropriate. 

An  alternate  approach  is  to  utilize  the  columns  of  the 
matrix  individually  to  discriminate  among  the  modal  components 
of  the  signal.  Since  all  components  of  a  given  modal  steering 
vector  share  a  common  value  of  these  phase  terms  may  be 
discarded  from  the  individual  elements  of  the  steering  vector. 
As  a  result  of  this  simplification,  the  modal  steering  vector 
is  completely  determined  by  the  (obtained  from  normal  mode 
modeling  of  the  immediate  environment)  and  the  desired  look 
direction  (determined  by  array  geometry) .  Having  established 
the  methodology,  and  replacing  the  phase  weights  with  steering 
delays,  the  output  of  the  modal  beamformer  is 

where  the  subscript  (m)  indicates  that  the  desired  output  is 
the  component  of  signal  representing  the  m^^  normal  mode. 
Expanding  this  form  to  explicitly  include  the  terms 
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representing  the  signal  and  noise  coinponents  of  the  acoustic 
field  yields 

r*i  n-l 

where  and  represent  the  amplitudes  of  signal  and  noise 
at  the  n^*^  hydrophone.  The  delay  and  weighting  of  the  noise 
received  at  the  hydrophones  has  no  undesirable  effects.  If 
the  noise  between  any  two  hydrophones  was  uncorrelated  prior 
to  the  operation,  then  it  will  remain  so  after  weighting. 
Again,  the  degree  to  which  the  ambient  noise  field  is 
uncorrelated  will  ultimately  reflect  upon  the  performance  of 
the  beamf ormer . 
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IV.  THE  HEARD  ISLAND  VERTICAL  ARRAY 


A.  ARRAY  CONSTRUCTION 

The  array  deployed  for  this  experiinent  consisted  of  32 
hydrophones  deployed  vertically,  each  having  a  nominal 
sensitivity  of  -170  dB  re  1  V/pPa.  The  element  spacing  was  45 
meters,  with  hydrophone  number  one  occupying  a  design  depth  of 
345  meters.  The  nominal  design  depth  for  hydrophone  number  32 
was  1740  meters. 

Instrumentation  on  the  array  consisted  of  two  sensors. 
The  upper  sensor  was  located  4.0  meters  above  hydrophone 
number  one.  It  recorded  ambient  pressure,  temperature,  tilt 
and  current  velocity.  The  lower  sensor  was  located  5.0  meters 
below  hydrophone  number  20.  It  recorded  tilt  (including 
direction),  pressure,  temperature  and  conductivity.  This 
sensor  appears  to  have  suffered  a  casualty  during  array 
deployment  which  rendered  its  tilt  data  suspect. 

Flotation  was  provided  by  one  main  syntactic  float  and  28 
syntactic  football  floats.  The  design  called  for  the  main 
float  to  reside  at  a  nominal  depth  of  23C  meters. 
Approximately  half  of  the  football  floats  were  submerged. 
This  arrangement  rendered  the  array  neutrally  buoyant,  thus 
providing  a  measure  of  isolation  from  surface  wave  effects. 
The  array  was  directly  tethered  to  the  R/V  Point  Sur  with 
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floatation  devices  on  the  surface  portion.  Each  float  was 
equipped  with  a  flasher,  a  flag  and  a  radar  reflector.  The 
float  nearest  the  ship  had  the  addition  of  a  radio  beacon. 
The  procedure  called  for  the  ship  to  remain  approximately  1200 
meters  from  the  position  immediately  above  the  array.  During 
the  experiment  the  R/V  Point  Sur  was  configured  to  remain 
quiet  and  dead  in  the  water,  except  as  required  to  keep  clear 
of  the  array.  Figure  (3.1)  illustrates  the  general 
configuration . 


Figure  3.1:  Heard  Island  Receiving  Array 
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B .  INSTRUMENT  DATA 


Ideally  the  receiving  array,  when  deployed,  is  oriented 
vertically.  However,  due  to  currents,  the  requirement  to 
maintain  control  of  the  research  vessel  and  other  factors 
beyond  the  control  of  the  research  team,  the  array  is  often 
tilted  from  the  vertical.  Given  the  length  of  the  array,  a 
1.0°  tilt  results  in  a  24  meter  horizontal  displacement 
between  hydrophones  one  and  32.  This  displacement  is  nearly 
one  (carrier  frequency)  wavelength.  Additionally,  the 
direction  in  which  the  array  tilts  affects  the  component  of 
horizontal  displacement  which  is  collinear  to  the  propagation 
direction.  Horizontal  displacements  perpendicular  to  the 
direction  of  propagation  are  acceptable  if  one  assumes  that 
the  wave  field  is  fairly  uniform  is  this  direction. 
Therefore,  the  array  steering  utilized  in  this  work  seeks  only 
to  correct  for  the  component  of  displacement  in  the  direction 
of  signal  propagation. 

The  speed  with  which  the  wave  field  propagates  is  required 
prior  to  the  calculation  of  steering  delays.  The  appropriate 
speed  to  use  in  this  case  is  the  group  speed  of  the  carrier 
frequency  for  the  mode  in  question,  defined  as 


^  d(ii 


(3.10) 


where  c^^  is  the  local  group  speed  for  the  m^^  normal  mode,  w 
is  the  radian  frequency  and  is  the  appropriate  eigenvalue. 
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The  amplitude  weight  applied  to  a  given  hydrophone  is  a 
function  of  the  mode  number  sought  and  hydrophone  depth.  The 
mode  shapes  {Z^(,z))  are  required  prior  to  beamforming.  The 
numerical  algorithm  used  to  obtain  the  eigenvalues  and 
eigenfunctions  for  this  work  was  developed  by  Chiu  and  Ehret 
[Ref.  5].  It  utilizes  finite  differencing  of  the  propagation 
vectors  for  a  given  frequency  and  sound  speed  profile  to 
calculate  the  modal  structure  of  the  local  environment.  The 
sound  speed  profiles  were  measured  during  the  conduct  of  the 
experiment  by  the  research  team  aboard  the  R/V  Point  Sur. 
Figure  3.2  is  the  sound  speed  profile  used  for  this  analysis. 
Figures  3.3  through  3.5  are  the  output  of  the  normal  mode 
model.  Amplitude  weights  for  the  array  elements  were 
calculated  from  these  data  sets. 
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Sound  Speed  Profile 


Sound  Speed  (meters/second) 

Figure  3.2:  Sound  speed  profile  at  receiving  array. 
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Depth  (meters) 


Mode  One 


Mode  Two 


Figure  3.3:  Normal  modes  one  and  two. 
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Depth  (meters) 


Mode  Three 


Mode  Four 
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Dopth  (meters) 


Mode  Five 


Mode  Six 
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C.  DATA  ACQUISITION  AND  PREPROCESSING 

Once  deployed,  the  receiving  array  recorded  numerous  data 
sets,  each  lasting  in  excess  of  60  minutes.  Each  recording 
comprises  32  channels  of  acoustic  data.  The  hydrophone  output 
was  ported  to  a  patch  panel  and  preamplifier  where  a  fixed 
gain  of  26  dB  was  applied.  Following  this,  the  signal  passed 
through  a  variable  gain  amplifier  which  automatically  set  the 
signal  gain  to  optimize  the  10  volt  dynamic  range  of  the 
analog  to  digital  converter.  This  process  had  the  duel  effect 
of  enhancing  the  precision  of  the  data  and  preventing 
saturation  of  the  A/D  converters.  Subsequent  normalization  of 
the  data  set  to  hydrophone  output  levels  was  accomplished  by 
an  implementation  of  a  program  written  by  Keith  Von  der  Heydt 
of  Woods  Hole  Oceanographic  Institution. 

Analog  to  digital  conversion  utilized  a  sampling  frequency 
of  228  Hz.  The  signal  was  passed  through  an  analog  bandpass 
filter  prior  to  conversion.  This  prevented  aliasing  of 
frequencies  above  Nyquist  into  the  data  set.  Low  frequency 
noise  components  were  also  attenuated.  The  cutoff  frequencies 
for  the  bandpass  filter  were  10  and  80  Hz.  Data  storage 
utilized  both  optical  disk  and  magnetic  tape.  Data 
acquisition  and  preprocessing  of  the  data  sets  used  in  this 
work  was  conducted  by  Miller,  Chiu,  Frogner  and  others. 
Figure  3.6  illustrates  the  equipment  alignment.  [Ref.  8] 
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Figure  3.6:  Data  acquisition  system 
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V.  RESULTS  AND  CONCLUSIONS 


The  desired  result  of  this  research  is  the  development  of 
a  software  package  capable  of  beamforming  a  tomography  signal. 
The  beamformer  should  discriminate  among  the  various  modal 
components  of  the  signal,  thus  enabling  one  to  treat  the  modes 
individually.  Due  to  the  inherent  nature  of  vertical  arrays, 
the  software  needs  to  incorporate  information  pertaining  to 
the  time  varying  array  tilt  and  individual  hydrophone  depths. 
To  this  end,  the  beamformer  utilizes  a  duty  cycle.  Each  such 
cycle  commences  with  an  assessment  of  array  geometry.  This 
information  is  used  to  determine  the  appropriate  steering 
delays  and  hydrophone  weighting  coefficients.  Once 
established,  the  beamformer  processes  60  seconds  of  acoustic 
data  before  repeating  the  procedure.  In  this  manner,  the 
software  package  implements  a  quasi-stable  virtual  array  in 
which  the  virtual  aperture  remains  fixed  in  space. 

The  results  presented  examine  two  performance  aspects. 
The  first  pertains  to  array  geometry  description.  Calculation 
of  the  steering  delays  and  amplitude  weights  is  detailed. 
Additionally,  the  acoustic  output  of  the  array  is  compared  to 
that  of  a  single  hydrophone  on  the  sound  channel  axis,  thus 
allowing  for  a  quantitative  evaluation  of  the  ultimate 
performance  of  the  beamformer. 
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A.  HYDROPHONE  DELAY  AND  WEIGHTING 


Inforroation  regarding  array  geometry  was  provided  by  two 
instrument  packages.  The  upper  sensor  was  located  4.0  meters 
above  hydrophone  number  one,  the  lower  sensor  5.0  meters  below 
hydrophone  number  20.  Figure  5.1  illustrates  the  instrument 
configuration  and  coordinate  system  used  in  determining  the 
array's  spatial  orientation. 


NORTH 


Figure  5.1:  Array  Instrumentation 
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1.  The  Upper  Sensor 

Information  provided  by  the  upper  instrument  includes 
the  component  of  current  velocity  in  both  the  north-south  and 
east-west  directions,  ambient  pressure,  temperature  and  tilt. 
The  direction  in  which  the  array  tilts  is  determined  by  the 
current  direction.  The  primary  assumptions  being  that  the 
current  acting  on  the  array  is  independent  of  depth  and  that 
the  direction  in  which  the  array  tilts  is  completely 
determined  by  the  current  direction. 

In  making  these  assumptions,  one  considers  the  primary 
source  of  current  measured  by  the  instrument  to  be  a  result  of 
ship  positioning  adjustments.  As  mentioned  previously,  the 
R/V  Point  Sur  was  required  to  periodically  move  the  ship  in 
order  to  remain  clear  of  the  array.  Such  movements  were 
intended  to  straighten  the  tether  between  the  array  and  the 
towing  rig,  thus  preventing  the  tether  from  becoming  fouled  in 
the  ship's  propellers.  These  adjustments  resulted  in  a 
tensioning  of  the  tether.  It  is  this  tensioning  that  is 
considered  to  have  resulted  in  array  movement  and  consequent 
water  flow  across  the  instrument  package.  Based  on  this 
assessment,  the  assumption  that  tilt  direction  was  determined 
by  the  flow  measured  at  the  upper  instrument  is  consistent 
with  the  conditions  of  the  experiment.  Higher  order  models 
could  be  derived  for  array  dynamics,  however  this  research  has 
utilized  the  first  order  model  of  a  linear  array. 
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The  next  variable  pertinent  to  array  geometry  is  the 
pressure  measured  by  the  upper  sensor.  From  this  it  is 
straight  forward  to  calculate  the  sensor’s  depth  using  the 
following  relationship, 

^  _  P  - 1  a  two  sphere 
P9 

where  P  is  the  pressure  (in  pascals)  measured  by  the 
instrument,  p  is  the  density  of  seawater  and  g  is  the  force  of 
gravity. 

Figures  5.2  and  5.3  present  the  data  received  from  the 
upper  instrument  package  during  two  of  the  acoustic  recording 
periods.  The  recording  beginning  00:05  (GMT)  January  27 
coincides  with  the  reception  of  a  m-seguence  of  length  2047. 
A  continuous  wave  signal  was  received  during  the  recording 
beginning  15:05  (GMT)  January  27.  They  are  representative  of 
the  most  and  least  favorable  conditions  encountered  during  the 
course  of  the  experiment. 
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Figure  5.2:  Array  geometry  for  00:05  (GMT)  January  27. 
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2.  The  Lower  Sensor 


The  lower  sensor  was  intended  to  provide  a  redundant 
source  of  array  positioning  data.  Specifically,  it  measured 
tilt  in  the  north-south  and  east-west  directions,  temperature, 
pressure  and  conductivity.  As  discussed  previously,  the  lower 
instrument  appears  to  have  suffered  a  casualty  which  rendered 
the  tilt  data  suspect. 

Figure  5.4  indicates  that  an  event  occurred  which 
caused  the  sensor  to  provide  data  which  is  not  consistent  with 
the  assumed  array  geometry.  One  scenario  explaining  the  event 
is  that  the  sensor  failed  shortly  after  entering  the  water. 
Another  possible  explanation  is  that  the  array  became 
entangled  in  itself  upon  descent. 

All  array  tilt  data  was  stored  for  subsequent  analysis 
ashore.  Therefore,  the  research  team  had  no  knowledge  of  the 
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3 .  Steering  Delays 


Time  delays  must  be  applied  to  each  of  the  hydrophones 
in  order  to  correct  for  array  tilt.  However,  only  the 
component  of  horizontal  displacement  collinear  to  the  signal 
propagation  direction  is  used  to  calculate  the  time  delay. 
This  requirement  mandates  apriori  knowledge  of  the  propagation 
direction.  Prior  work  has  determined  that  the  expected 
direction  from  which  the  signals  would  arrive  is  217.0'^  true 
[Ref.  9]. 

In  addition  to  array  geometry  and  propagation 
direction,  one  requires  a  nominal  propagation  speed  to  convert 
the  horizontal  displacements  to  time  delays.  The  appropriate 
speed  to  use  is  the  mode  group  speed,  as  previously  discussed. 
Group  speeds  for  this  study  were  obtained  from  eigenvalues 
calculated  by  the  normal  mode  model  courtesy  of  Chiu  and  Ehret 
[Ref.  5].  The  mode  group  speed  of  the  carrier  frequency  was 
used  to  calculate  time  delays.  Figure  5.5  presents  the  mode 
group  speeds  as  a  function  of  frequency  for  the  lowest  five 
modes . 
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Figure  5 <,5:  Mode  group  speeds. 
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4.  Hydrophone  Amplitude  Weighting 

Application  of  hydrophone  weights  utilizes  the  array 
geometry  to  determine  the  depth  of  each  array  element.  The 
value  of  the  depth  dependent  function  is  then  assigned  to 
the  appropriate  hydrophone  output.  This  manner  of  weighting 
enhances  the  array  sensitivity  to  the  m*^  normal  mode. 

An  additional  consideration  encountered  in  hydrophone 
weighting  is  fault  tolerance ,  Vertical  acoustic  arrays 
operate  in  a  hostile  environment,  thus  individual  hydrophones 
may  fail.  The  amplitude  weights  applied  to  the  array  must 
cope  with  the  failure  of  individual  elements.  To  this  end, 
the  beamformer  applies  a  weight  of  zero  to  elements  which  have 
been  evaluated  as  unreliable. 

Log  books  kept  by  the  research  team  aboard  the  R/V 
Point  Sur  indicating  suspect  hydrophones  were  verified  by  post 
processing  statistical  analysis  of  the  individual  elements. 
Several  hydrophones  were  found  to  have  failed  during  the 
course  of  the  experiment.  Specifically  elements  18  through  28 
were  evaluated  as  unreliable.  As  a  result  of  this 
determination,  only  hydrophones  one  through  17  were  used  to 
beamform  the  target  signals. 

B.  ACOUSTIC  PERFORMANCE 

The  preliminary  data  analysis  seeks  to  quantify  the 
beamformer 's  performance  with  respect  to  a  single  hydrophone 
located  on  the  sound  channel  axis.  The  power  spectrum 
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calculated  from  the  output  of  hydrophone  five  was  compared  to 
the  spectra  calculated  from  the  beamformed  output  of  channels 
one  through  17.  The  array  gain  is  defined  as 

where  SNR^  is  the  signal  to  noise  ratio  of  the  beamformed 
output  and  SNR^  is  the  signal  to  noise  ratio  of  the  n'^ 
hydrophone  [Ref.  10].  The  SNR  for  each  case  was  calculated  by 
integration  of  the  power  spectra  over  a  6.0  Hz  bandwidth 
centered  on  the  continuous  wave  frequency. 

The  array  gain  was  calculated  for  numerous  power  spectra / 
each  representing  60  seconds  of  acoustic  data.  Additionally, 
the  gain  was  calculated  for  the  lowest  five  normal  modes 
present.  Values  for  array  gain  range  from  a  low  of  2.5  dB  to 
a  high  of  8.5  dB.  The  nominal  array  gain  is  6  dB. 

Reasons  for  the  range  in  values  include  temporal 
variability  in  the  modal  components  of  the  signal,  errors 
introduced  in  the  steering  delays  and  errors  in  amplitude 
weighting.  Time  delay  and  amplitude  weight  errors  are  the 
direct  result  of  assumptions  made  regarding  array  geometry. 
The  primary  source  of  geometric  error  is  considered  to  involve 
the  assumption  that  tilt  direction  is  completely  determined  by 
current  direction.  Array  tilt  bearing  may  not  be  aligned  well 
with  current  direction  during  periods  when  the  direction  of 
flow  across  the  sensor  is  varying.  Therefore,  during  periods 
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(V'/Hz) 


of  variable  current  flow,  marginal  performance  can  be  expected 
of  the  beamformer. 


Figures  5.6  through  5.11  present  representative  power 
spectra  of  single  channel  data  and  beamformed  data  for  modes 
one  through  five  over  the  entire  range  of  frequencies 
analyzed . 
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Figure  5.6:  Hydrophone  five  power  spectrum. 
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Figure  5.8:  Mode  two  power  spectrum. 
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Mode  three  power  spectrum. 


O  20  40  60  80  100  120 

Frequency  (Hz) 

Figure  5.10:  Mode  four  power  spectrum. 
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Figure  5.11:  Mode  five  power  spectrum. 


As  shown  in  the  previous  figures,  the  noise  field  at  the 
beamformer  is  not  incoherent.  This  fact,  coupled  with 
imperfect  geometric  array  description,  detracts  from  the 
calculated  array  gain. 

One  feature  worth  noting  is  the  decreasing  strength  of  the 
undesirable  tonal  components  with  increasing  mode  number, 
specifically  those  in  the  40  to  50  Hz  range.  Presumably, 
these  components  were  radiated  by  the  R/V  Point  Sur.  One 
hypothesis  is  that  the  strength  of  these  signals  in  mode  one 
is  the  result  of  grating  lobes  in  the  array's  beam  pattern  for 
this  mode.  The  amplitude  weights  applied  to  the  hydrophones 
are  similar  to  those  for  a  plane  wave  beamformer,  in  that 
there  is  no  phase  shift  applied  along  the  array  in  the  form  of 
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coefficients  of  differing  sign.  Therefore,  spatial  aliasing 
is  possible  for  frequencies  with  wavelengths  less  than  twice 
the  hydrophone  spacing.  For  the  Heard  Island  array,  this 
includes  frequencies  ranging  from  17  Hz  to  the  bandpass  cutoff 
frequency  of  80  Hz. 

As  the  mode  number  sought  increases,  these  frequency 
components  diminish  in  strength.  This  indicates  that  the 
rejection  of  coherent  noise  of  local  origin  improves  as  the 
number  of  180°  phase  changes  applied  to  the  hydrophone  output 
increases . 

Also  shown  in  these  figures  is  the  presence  of  distant 
sources  of  coherent  noise.  These  other  tonal  components  are 
seen  to  vary  both  as  a  function  of  mode  number  and  time.  This 
observed  behavior  is  consistent  with  modal  propagation  in  a 
range  dependent  environment.  Merchant  traffic  in  shallow 
water  and  offshore  drilling  activity  are  two  possible 
candidates  for  sources  of  coherent  noise  which  becomes  coupled 
into  the  sound  channel. 

Figures  6.12  through  5.17  present  the  same  power  spectra 
in  the  immediate  vicinity  of  the  carrier  frequency  (57  Hz). 
Again,  the  beamformed  output  indicates  successful  attenuation 
of  undesirable  components  of  a  coherent  noise  field, 
presumably  of  local  origin. 
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Power  Spectral  Density  (VV  Hz)  Power  ^pecM  Density  (V^/ Hz) 


Frequency  (Hz) 

lire  5.12:  Hydrophone  five  power  spectrum. 


Frequency  (Hz) 

Figure  5.13:  Mode  one  power  spectrum. 
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Mode  Two  Power  Spectrum 


Mode  Four  Power  Spectrum 


Frequency  (Hz) 

Figure  5.16:  Mode  four  power  spectrum. 
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C.  CONCLUSIONS 


Modal  beamforming  is  a  valuable  technique  for  use  in  long 
range  tomography  experiments.  The  signals  presented  ir  this 
work  traversed  the  longest  transmission  path  for  which 
reception  was  attempted  during  the  course  of  the  Heard  Island 
Experiment.  Although  no  travel  time  estimates  have  been  made, 
the  successful  detection  of  the  signals  indicates  that  large 
scale  tomography  experiments  are  feasible. 

The  primary  concern  of  this  work  was  the  development  of 
a  system  capable  of  detecting  the  acoustic  signals  emanating 
from  the  vicinity  of  Heard  Island.  This  has  been 
accomplished.  Additionally,  the  data  describing  the  complex 
dynamics  of  the  near  vertical  acoustic  array  has  implications 
lor  the  design  and  construction  of  vertical  arrays  for  use  in 
future  tomography  experiments.  The  most  troublesome  aspect  of 
this  work  has  been  the  accurate  determination  of  array 
geometry,  specifically  tilt  direction.  Describing  tilt 
direction  in  terms  of  measured  current  at  one  point  on  the 
array  is  considered  inadequate. 

D .  RECOMMENDAT I ONS 

Improvements  to  this  implementation  include  physical 
aspects  of  the  array  itself  and  performance  of  the  numerical 
algorithm.  Large  vertical  arrays  are  particularly  sensitive 
to  errors  in  estimated  geometry.  In  this  experiment,  an  error 
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of  0.25°  in  tilt  (assuming  perfect  knowledge  of  tilt 
direction)  creates  a  phase  error  of  45  degrees  between  the 
most  distant  elements.  This  error  is  more  destructive  tc  che 
higher  modes  owing  to  their  greater  spatial  extent.  The 
problem  is  aggravated  by  questionable  tilt  direction  data. 

Construction  of  future  vertical  arrays  should  incorporate 
multiple  instrument  packages  which  directly  measure  tilt  and 
direction.  As  this  was  the  case  with  the  lower  instrument  on 
the  Heard  Island  array,  it  is  most  unfortunate  that  it  failed. 
Additionally,  an  investigation  into  the  utility  of  higher 
order  models  describing  array  shape  would  be  illuminating. 
Incorporating  a  quadratic  model  for  array  shape  would  not 
impose  a  significant  computational  load  on  the  beamformer  and 
may  improve  the  estimated  position  of  individual  hydrophones. 

The  current  implementation  of  this  beamfoimer  utilizes 
third  order  polynomial  interpolation  during  the  application  of 
time  delays  to  individual  hydrophone  outputs.  As  such  it 
represents  the  greatest  computational  load  in  the  program. 
The  interpolator  is  based  on  Neville's  algorithm  [Ref.  11]. 
No  attempt  was  made  to  minimize  execution  time.  As  a  result, 
the  software  requires  approximately  four  hours  to  process  one 
hour  of  acoustic  data.  A  more  efficient  interpolator  may 
reduce  execution  time  sufficiently  to  permit  this  technique  to 
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be  used  for  real  time  modal  beamforming.  To  this  end,  the 
software  is  assembled  from  modular  components,  thus  allowing 
for  improvements  on  a  function  by  function  basis. 
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APPENDIX  A 


The  following  programs  were  developed  during  the  course  of 
this  research.  They  include  the  modal  beamformer  and 
associated  utility  programs.  These  programs  may  be  obtained 
by  contacting  James  H.  Miller  at  his  address  in  the  initial 
distribution  statement. 

A.  THE  MODAL  BEAMFORMER 

1.  Operational  Considerations 

The  modal  beamformer  was  designed  with  portability  a 
prime  consideration.  As  a  result,  there  are  numerous  program 
parameters  which  must  be  set  by  the  end  user.  Among  these  are 
the  characteristics  of  the  physical  array  and  data  acquisition 
system.  All  such  user  selectable  parameters  are  contained 
solely  within  the  program  definitions.  The  parameters  appear 
below  in  the  same  order  in  which  they  appear  in  the  program. 

•  UNIX  VERSION:  selectively  compiles  either  unix  or  ansi 
compatible  code  blocks. 

•  ASCII:  selects  output  in  either  ascii  or  binary  format. 
Ascii  is  useful  for  data  export,  binary  saves  storage 
space . 

•  SIGNAL:  permits  user  to  enable  or  disable  the  beamformer 's 
time  series  output. 

•  SPECTkUM:  permits  the  user  to  enable  or  disable  the 
beamformer 's  power  spectrum  output. 
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•  LOWER_SENSOR:  directs  the  program  to  load  or  ignore  lower 
tilt  sensor  data.  This  option  should  remain  OFF  as  there 
is  no  reliable  lower  sensor  data. 

•  VALIDATE:  enables  program  validation.  If  selected  the 
program  will  create  output  files  which  list  the  array 
geometry,  steering  delays,  hydrophone  weights  and 
calculated  mode  group  speed.  This  feature  is  useful  for 
program  debugging. 

•  ERROR_ESTIMATE :  if  selected,  this  feature  causes  the 
polynomial  interpolator  to  store  the  upper  bound  on  the 
mean  squared  error  in  the  steering  delays.  This  feature 
doubles  the  output  of  the  beamformer. 

•  ON:  logical  switch  for  program  control. 

•  OFF:  logical  switch  for  program  control. 

•  INTERPOLATE:  selects  interpolator.  This  feature  permits 
the  installation  of  an  improved  interpolator  without  the 
requirement  to  search  for  and  replace  each  call  to  the 
function . 

•  ORDER:  indicates  the  degree  of  the  polynomial  used  in  the 
polynomial  interpolator. 

•  STEP:  indicates  the  number  of  points  on  either  side  of  a 
sequence  to  be  taken  when  estimating  derivatives. 

•  TINY:  prevents  division  by  zero  in  floating  point 
operations . 

•  PI:  used  for  trigonometric  recursions. 

•  RADIAN:  used  for  degree  to  radian  conversions. 

•  OFFSET:  indicates  the  distance  (in  meters)  between  the 
number  one  hydrophone  and  the  upper  tilt  sensor. 

•  DELTA_R:  indicates  array  element  spacing  (in  meters). 

•  CTD_OFFSET:  indicates  the  difference  in  depth  between  the 
sound  speed  profile’s  first  data  point  and  the  depth 
increment  used  in  the  profile  (in  meters). 

•  SSP_LENGTH:  indicates  the  maximum  number  of  data  points 
that  an  eigenfunction  may  contain.  The  only  restriction 
applicable  to  this  parameter  is  available  memory. 
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•  EIGVAL_LENGTH :  indicates  the  itiaxiraurc  number  of 

eigenvalues  that  the  program  will  accept.  The  only 
restriction  applicable  to  this  parameter  is  available 
memory . 

•  LOOK_DIRECTION:  the  compass  direction  (in  degrees  true) 
from  which  the  signal  arrives. 

•  TILT_BUFFER:  the  buffer  length  allocated  to  tilt  data. 
The  only  restriction  applicable  to  this  parameter  is 
available  memory. 

•  BUFFER_TIME:  the  length  of  time  (in  seconds)  represented 
by  the  acoustic  data  input  buffer.  This  period  represents 
the  output  buffer  length  plus  one  second  on  either  end  to 
permit  steering  delays  and/or  advances.  The  array  is 
currently  capable  of  end  firinc^a  1600  meter  acoustic 
array.  This  parameter  must  be  aninteger.  Additionally, 
60  must  be  an  integer  multiple  of  (BUFFER_TIME-2 ) .  The 
array  has  a  duty  cycle  of  one  minute  between  design 
changes.  There  must  be  an  integer  number  of  output 
buffers  per  duty  cycle. 

•  F_SAMPLE:  sampling  rate  of  the  data  acquisition  system. 
This  value  must  be  an  integer. 

•  CHANNELS:  the  number  of  hydrophones  on  the  array. 

•  F_CARRIEL:  the  carrier  frequency  of  the  target  signal. 
The  group  speed  used  in  calculating  steering  delays  is 
based  on  this  frequency.  This  must  be  entered  in  floating 
point . 

•  FFT_LENGTH:  the  length  of  the  FFT  used  if  frequency 

domain  output  is  selected.  This  value  must  be  a  power  of 
two  which  is  less  than  or  equal  to  the  number  of  samples 
in  a  single  channel  of  input  data. 

•  SWAP:  a  macro  used  in  calculating  a  FFT. 


The  following  program  functions  are  adaptations  of 
those  found  in  Reference  11: 


•  polint: 

•  realft: 

•  fourl: 


performs  the  polynomial 
calculates  Fast  Fourier 
calculates  Fast  Fourier 


interpolation . 
Transforms . 
Transforms 
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•  vector:  allocates  memory  for  one  dimensional  arrays. 

•  matrix:  allocates  memory  for  two  dimensional  arrays. 

•  free_matrix:  deallocates  memory  from  two  dimensional 

arrays . 

•  ExitOnError:  provides  abnormal  program  termination. 


2.  Beamformer  Source  Code 

/*★** 

*  PROGRAM:  BEAMFORMER  \/sn  1.0 

*  WRITTEN  BY;  Steven  Crocker 

*  LAST  UPDATE:  October  9,  1991 

* 

*  This  program  takes  input  from  various  data  files  and  the  user.  It 

*  outputs  a  data  file.  The  inputs  are  a  number  of  channels  of  digital 

*  acoustic  data,  and  information  regarding  the  physical  characteristics 

*  and  geometry  of  the  receiving  array.  Additionally,  environmental  data 


*  in  the  form  of  normal  mode  eigenfunctions  and  eigenvalues  at  the 

*  receiving  array  are  required  to  operate  this  beamformer.  The  output 

*  is  a  single  channel  of  acoustic  data. 

* 

(Kdefine  UNIX  VERSION  /*  either  ANSI  or  UNIX  */ 

^define  ASCII  ON  /*  select  output  mode  ON  or  OFF  */ 

((define  BINARY  OFF  /*  select  output  mode  ON  or  OFF  */ 

((define  SIGNAL  ON  /*  either  ON  or  OFF  */ 

((define  SPECTRUM  ON  /*  either  ON  or  OFF  */ 

((define  LOWER  SENSOR  OFF  /*  either  ON  or  OFF  */ 

((define  VALIDATE  OFF  /*  either  ON  or  OFF  */ 

((define  ERROR_ESTIMATE  OFF  /*  either  ON  or  OFF  */ 

((define  ON  1  /*  logical  "switch"  */ 

((define  OFF  0  /*  logical  "switch"  */ 

((define  INTERPOLATE  polint  /*  polint  */ 

((define  ORDER  3  /*  order  of  interpolator  (odd)*/ 

((define  STEP  1  /*  number  of  steps  for  derivatives  */ 

((define  TINY  1.0e-25  /*  prevents  division  by  zero  .  */ 

((define  PI  3.14159265359  /*  for  freq  to  omega  conversions  */ 

((define  RADIAN  57.2957795131  /*  for  degree  to  radian  conversions  */ 

((define  OFFSET  4.0  /*  dist  btwn  upper  sensor  and  phone  ((1  */ 

#define  DELTA_R  45.0  /*  array  element  spacing  */ 

#define  CTD_0FFSET  0.5  /*  diff  btwn  ctd  depth  inc  S  1st  depth  */ 

((define  SSP_LENGTH  2500  /*  max  number  of  pts  in  eigunfunction  */ 

((define  EIGVAL_LEN6TH  230  /*  max  number  of  eigenvalues  */ 

((define  L00K_DIRECTI0N  217.0  /*  direction  from  which  signal  arrives  */ 

((define  TILT_BUFFER  120  /*  max  length  of  tilt  data  vectors  */ 

((define  BUFFER_TIME  12  /*  input  buffer  length  in  seconds  (int)*/ 

((define  F_SAMPLE  228  /*  sampling  frequency  (int)*/ 

#define  CHANNELS  32  /*  number  of  channels  processed  */ 

((define  F  CARRIER  57.0  /*  carrier  frequency  */ 

((define  FFT_LENGTH  2048  /*  radix  2  <=  (BUFFER_T1ME-2)*F_SAMPLE  */ 

•define  SUAP(a,b)  tempr=(a); (a)=(b); (b)=tempr 


((include<stdio.h> 

#include<malloc.h> 

((include<math.h> 

((if  defined  (  ANSI  ) 
((include<f  loat .  h> 
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#include<stdlib.h> 
int  getInputCvoid); 
int  putOutput(void); 

int  processTi It (float  **x,  float  **y,  float  **2); 

int  processHodesC float  **z,  float  *weight,  float  *ptrC); 

int  dydxCfloat  *x,  float  *y,  float  *ddx,  int  points); 

int  polintCfloat  *xa,  float  *ya,  int  n,  float  x,  float  *y,  float  *dy); 

int  realftCfloat  *data,  int  n,  int  isign); 

int  fourl (float  *data,  int  nn,  int  isign); 

float  *vector(int  Length); 

float  **iiatrix(int  row,  int  col); 

int  f ree_«atrix(f loat  **in,  int  row); 

void  ExitOnError(ehar  error_txt[]); 

<Kelif  defined  (  UNIX  ) 
int  vectorO, 

•atrixO; 

get Input (),  • 

putOutputO, 
processTi  ItO, 
processNodesO, 
dydx ( ) , 
polintO, 
realftO, 
fourl {), 
free_inatrix(), 

ExitOnErrorO; 

#endif 

/*  Global  Variable  Declarations  */ 
float  **inSOUND,  *outS0UND,  *MeanSqError,  Max=0.0,  f1in=1.0e25; 
int  lastMinute,  Hinute=1,  if irstBuffer=1; 

mainO 

{ 

FILE  *fpV1; 

int  i,  k,  n,  *shift; 

float  **x,  **y,  **z,  *indx,  *samples,  arg,  ans,  err,  *delay, 

♦delta,  *weight,  *pwrSpectrum,  Cgroup; 


/♦  Memory  Allocation  and  Tilt  Data  Processing  */ 
x=(f loat**)matrix(CHANNELS,  TILT_BUFFER); 
y=(f loat**)matrix(CHANNELS,  TILT  BUFFER); 

2=(floate*)matrix( CHANNELS,  TILt'buFFER); 
processTi lt(x,  y,  2); 

if  (VALIDATE==ON) 

{ 

printf ("\nDumping  array  geometry  to  array.datXn"); 
f pVI =f open ("array. dat " , "wt " ) ; 

if (fpV1==NULL)  ExitOnErrorC'Error  opening  validation  file"); 

fprintf (fpVl,"Channel\t\t  X\t\t  Y\t\t  Z\n"); 

for  (i=1;i<=lastMinute,i++) 

{ 

fprintf(fpVl,  "MINUTE:  Xi\n",  i); 
for  (j=1; j<=CHANNELS; j++) 

fprintf(fpVl,"%i\t  /ifXt  Xf\t  %f \n", j,xtj]Ci3,y[j]Ci],2[j]Ci]); 
>  /♦  for  ♦/ 
f close (fpVl ); 

>  /*  if  */ 

/*  Memory  Allocation*/ 
weighted  loat*)vector(  CHANNELS); 
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1ndx=(float*)vector(0RDER+1); 
delBy= (float*) vectorC CHANNELS); 
delta=(f loat*)vector(CHANNELS); 
outSOUND=(float*)vector((BUFFER  TIHE-1)*F_SAHPLE); 
inSOUND=(f Loat**)aatrix(CHANNELS,  BUFFER  T1ME*F_SAHPLE); 
if  ((shift=(int*)«iaUoc((CHANNELS+1)*sizeof(int)))==NULL) 
ExitOnErrorCMenory  allocation  failure  for  shiftC].”); 
for  (i=1;  i<=ORDER+1;  i++)  indxCiXf loat)i; 
if  (ERROR  ESTIMATE==ON) 

HeanSqError=(float*)veetor((BUFFER_TIHE-1)*F_SAMPL£); 
if  (SPECTRUM==ON)  pwrSpeetruB=(f loat*)vector(FFT_LENGTH/2); 
whi le(Minute<=lastMinute) 

{ 


/*  Node  Data  Processing  */ 
processNodes(z,  weight,  &Cgroup); 

/*  Calculate  delays,  shifts,  etc  */ 
for  (i=1;  i<=CHANNELS;  i++) 

< 

delay[i]=xCi][Hinute]/Cgroup;  /*  Time  delay  */ 

shift[i]=(int)(deLay[i]*(f loat)F_SAMPLE);  /*  #  of  samples  */ 

/*  fraction  of  1  sample  */ 
deltaCi]=delay[i]*(f loat)F_SAMPLE-(f loat)shift[i]; 

>  /*  for  */ 

if  (VALIDATE==ON) 

{ 

printf ("Dumping  element  delays  to  delay.dat\n"); 

if  (f irstBuffer) 

C 

if  ((fpV1=fopen("delay.dat",  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file."); 

>  /*  if  */ 

else 

{ 

if  ((fpVl*fopen("deLay .dat",  "at"))  *=  NULL) 

ExitOnErrorC'Error  opening  validation  file."); 

>  /*  else  */ 

fprintf(fpVl,  "HINUTE;  /iiVn",  Minute); 
fprintf (fpVl , 

"ChannelXt  delay\t\t  int  shiftXt  fraction  of  1  shiftVn"); 

for  (i=1;i<=CHANNELS;i++) 

{ 

fprintf (fpVI,  "r.i\t  ;;e\t  Xi\t:;  e\p", 
i,delay[i],shiftCi],deltaCi3); 

>  /*  for  */ 

fclose(fpV1 ); 

printf ("Dumping  phone  weights  and  grp  speed  to  modal .dat\n"); 

if (f irstBuffer) 

{ 

if  ((fpV1=fopen("modal.dat",  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file.Xn"); 
fprintf(fpVl, "Group  speed  for  F_CARRIER  is:  Xg\n\n\n",Cgroup); 

>  /*  if  */ 

else 

{ 

if  ((fpV1=fopen("modal.dat","at"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file.Vn"); 
fprintf  CfpVl ,  "\nHydrophone  weights  for  minute  5ii\n", Minute); 

>  /*  else  if  */ 
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fprintf (fpVl,  "Channel  \tUeight\n"); 

for  (i=1;i<=CHANNELS;i++) 

fprintf (fpVI,  "Xi\tXe\n",i,weight[i3); 

fclose(fpVI); 

>  /*  if  */ 

if  (SPECTRUH==ON) 

for  (k=1;k<=FFT_LEN6TH;k++)  pwrSpectruiBCk3=0.0; 

for  (n=1;  n<=60/(BUFFER  TIHE-2);  n++) 

{ 

get  Input  O; 

if (f irstBuffer)  /*  Produce  first  output  buffer  */ 

< 

for  (i=1;  i<=F  SAMPLE;  i++) 

{ 

out SOUNDCi 3=0.0; 

Tf  (ERROR  ESTIMATE==ON)  MeanSqErrorCi 3=0.0; 

}  /*  for  */  " 

for  (i=F_SAMPLE+1;  i<=(BUFFER_TIME-1)*F  SAMPLE;  i++) 

< 

outSOUNOCi 3=0.0; 

if  <ERR0R_EST1MATE==0N)  MeanSqErrorCi3=0.0; 

for  (j=1;  i<=CHANNELS;  j++) 

{ 

if  (delayC j3>=0.0) 

{ 

arg=(float)(ORDER+1)/2.0+(1-deUaCj3); 
samples  =  &inSOUNOCj3[i-shiftCj3-(ORDER*’1  )/2]; 

)  /*  if  */ 

else  if  (delayCj3<0.0) 

arg=(f loat)(0RDER♦1)/^.0+deltaCj3; 

samples  =  &inSOUNOCj3Ci-shiftCj3-(ORDER+1)/2+13; 

}  /*  else  if  */ 

INTERPOLATEdndx,  samples,  OROER+1,  arg,  Rans,  Serr); 
outS0UNDCi3=outS0UNDCi]+ans; 
if  (ERR0R_EST1MATE==0N) 

MeanSqEr ror [ i 3=MeanSqE  r ror[ i 3+err*err; 

>  /*  for  */ 

if (fabs (out SOUNDCi 3 )>Max)  Max=fabs(outS0UND[i3); 
i f ( f abs (out SOUND [ i 3 )<M i n )  Mi n=f abs(outSOUND[ i 3 ) ; 
if  (ERR0R_ESTIMATE==0N) 

MeanSqErrorCi ]=MeanSqErrorCi3/(f loat)CHANNELS; 

>  /*  for  */ 

)  /*  If  */ 

else  /A  Produce  subsequent  output  buffers  */ 

for  (i=F_SAMPLE+1;  i<=(BUFFER  TIME-1)*F  SAMPLE;  i++) 

C 

outSOUNDCi-F_SAMPLE]=0.0; 

if  (ERR0R_ESTIMATE==0N)  MeanSqErrorCi-F_SAMPLE3=0.0; 

for  (j=1;  j<=CHANNELS;  j++) 

{ 

if  (delayC j3>=0.0) 

t 

arg=(float)(ORDER+1)/2.0+(1-deltaCj3); 
samples  =  SinSOUNDC jlCi-shiftCj3-(0RDER+1 )/23; 

>  /*  if  */ 

else  if  (delayCj3<0.0) 

arg=(f lost ) (ORDER+1 )/2.0*deltaCj3; 

samples  =  8inS0UNDCi3Ci-shiftCj3-(0RDER+1)/2+13; 
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}  /*  else  if  */ 

INTERPOLATECindx,  samples,  OROER+1,  arg.  Sans,  Serr); 
ootSOUNOCi-F  SAt1PLE]=outS0UNDCi-F  SAHPLE]+ans; 
if  (ERROR_ESTIHATE==ON) 

HeanSqErrorCi-F  SA?1PLE]=HeanSqEprorCi-F  SAHPLE3+err*err; 

}  /*  for  */  ~ 

iflfabsCoutSOUNDCi-F  SANPLE])>nax) 

«ax=f abs (out SOONOC i -F_SAHPLE] ) ; 

if(fabs(outSOUND[i-F  SAHPLE])<nin) 

«in=fabs(ootSOONDri-F  SAMPLE3); 
if  (ERROR_ESTIMATE==ON)  ~ 

MeanSqirrorCi-F  SAHPLE]=HearSaErrorCi-F  SAMPLED/ 

( float /CHANNELS; 

}  /*  for  */ 

>  /*  else  */ 

if  (SIGNAL==ON)  putOutputO; 

if  (SPECTRUH==ON) 

{ 

windowCoutSOUNO,  FFT_LEN6TH); 
realftCoutSOUND,  FFT~LENGTH/2,  1); 
for  (k=0;k<FFT  LENGTH;k  +=  2) 

( 

purSpectruinCk/2+1]=pwrSpectrumCk/2+13+ 

outS0UNDCk]*outS0UND[k]+outS0UN0Ck+1D*outS0UND[k+13; 

}  /*  for  */ 

}  /*  if  */ 

f irstBuffer=0;  /*  set  firstBuffer  to  false  */ 

}  /*  for  */ 

if  (SPECTRUM==ON) 

{ 

for  (k=1;  k<=FfT_LENGTH;  k++) 

pwrSpectrumCk3=purSpeetrumCkD*(f loat)(BUFFER_TIME-2)/60.0; 

dufflpSpectrum(pwrSpectrum); 

}  /*  if  */ 

printf ("\t;'.i  minutes  of  input  data  processed. \n".  Minute); 
Minute*+;  /*  increment  minute  counter  */ 

>  /*  while  */ 

printf ("EXECUTION  COMPLETE;  End  of  tilt  data  encountered\n"); 
printf ("Maximum  magnitude  encountered  was;  %e\n",Max); 
printf ("Minumum  magnitude  encountered  was:  Xe\n",Min); 
exit(O); 


FUNCTION:  getInputO 

This  function  handles  all  acoustic  input.  It  also  provides  one  of 
two  normal  process  terminations  available  in  the  program.  (The  other 
is  located  in  mainO  . ) 


Arguments; 

Return  value; 

Functions  called: 
Definitions  called: 

Global  variables  called: 


none 

0 

vector () 
ANSI 

F  SAMPLE 

bUffer_time 

inSOUNOCDCD 

firstBuffer 


ExitOnErrorO 

UNIX 

CHANNELS 


Min 

Max 
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*  Significant  aemory  allocation:  diskBufferC3 

* 

*****/ 

iKif  defined  (  ANSI  ) 
int  getlnput(void) 
fDelif  defined  (  UNIX  ) 
get Input O 
(Kendif 
{ 

int  i,  j,  buffer,  items; 
float  AdiskBuffer; 
char  f i leNameEBO]; 
static  FILE  efpStatic; 

if  (firstBuffer) 

{ 

printf ("Enter  file  name  for  input  acoustic  data:  "); 

if((scanf("Xs",  file;(ame))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fpStatic=fopen(f i leName,  “rb"))  ==  NULL) 

ExitOnErrorC'Error  opening  INPUT  ACOUSTIC  data  file."); 
>  /*  if  */ 


/*  Hemory  Allocation  */ 

if (firstBuffer)  buf f er=BUFFER  TIME*F  SAMPLE*CHANNELS; 
else  buffer=(BUFFER_TIME-2)*F~SAMPLE*CHANNELS; 
diskBuf fer=( float*) vector (buffer); 


items=fread((ehar*)(diskBuffer+1),  sizeof (float),  buffer,  fpStatic); 
if  (items==buffer)  /*continue*/; 
else  if (ferror(fpSt8tic)  !=  0) 

ExitOnErrorC'Error  encountered  while  reading  input  acoustic  data") 
else  if (feof (fpStatic)  !=  0) 

{ 


printf (" 

printf ("\t  End  of  File  reached;  EXECUTION  CO«PLETE\n"); 
printf ("\t\t%i  minutes  of  data  processed\n",Hinute-1); 
printf ("\t\t‘/.i  bytes  of  data  discardedXn",  items*sizeof (float)); 
printf ("\t\tMaximum  magnitude  encountered  was:  %e\n",Hax); 
printf ("\t\tMinimum  magnitude  encountered  was;  Xe\n",Min); 
printf ("\t  End  of  File  reached:  EXECUTION  COMPLETE\n"); 

fclose(fpStatic); 

exit(O); 

>  /*  else  if  */ 

else  ExitOnError( "Unknown  error  handling  acoustic  input  file."); 


for  (i=1;  i<=BUFFER_TIME*F  SAMPLE;  i++) 

< 

for  (j=1;  j<=CHANNELS;  j++) 

if  (firstBuffer) 

{ 

inSOUNDCjICi]  =  diskBufferCCHANNELS*(i-1)+jD; 

}  /*  if  */ 

else 

{ 

if  (i<=2*f_SAfiPLE) 

inSOUND[j][i]=inSOUNOCi][i+(BUFFER_TIHE-2)*F_SAMPLE]; 

else 

inSOUNO[j][i]=diskBuffer[CHANNELS*(i-2*F  SAMPLE-I )+j]; 
>  /*  else  */ 

>  /*  for  */ 

>  /♦  for  */ 
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/*  Deallocate  Memory  */ 
f  reel  (ehar*)dislcBuf  f  er); 
returnl  0  ); 


/**♦**  end  get  Input  **★**/ 

/AAtAAAAAAAAAAAAAAAAAAAAA/ 

/*★*** 

*  FUNCTION:  putOutputO 

* 

This  function  handles  all  acoustic  output.  Additionally,  it  outputs 
the  estimated  mean  squared  error  from  the  interpolators  (if  enabled). 


Arguments: 

Return  value: 

Functions  called: 
Definitions  called: 

Global  variables  called: 


none 

0 

ExitOnErrorO 

ANSI 

F_SAMPLE 

outSOUNDC] 
f irstBuffer 


UNIX 

BUFFER_TIME 

MeanSqErrorC] 


Significant  memory  allocation:  none 


«if  defined  (  ANSI  ) 
int  putOutput(void) 

«elif  defined  (  UNIX  ) 

putOutputO 

Aendif 

< 

int  i,  cut; 

char  f i leNameC12],  modeC2]; 
static  FILE  AfpOutSound,  *fpMSE; 

if (f irstBuffer) 

{ 

if  (ASCII==ON) 

modeC0]= ' w ; 
modeCI ]- ' t  ; 

>  /*  if  */ 

else  if  (BINARY==ON) 

mode[0]='w' ; 
mode[1]=  'b'  ; 

>  /*  else  if  */ 

cut=1 ; 

printf ("Enter  file  name  for  output  data:  "); 

if ((scanf ("as",  f i leName) )==E0F) 

ExitOnErrorCFatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fpOutSound=fopen(fi leName,  mode))  ==  NULL) 
ExitOnErrorC'Error  opening  OUTPUT  data  file."); 

if  (ERROR_ESTIHATE==ON) 

printf ("Opening  file  error. dat\n\n"); 

if  ((fpnSE=fopen("error.dat",  mode))  ==  NULL) 
ExitOnErrorC'Error  opening  error.dat"); 

>  /*  if  */ 
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>  /*  if  */ 
else  cut=2; 


if  (ERROR  EST1HATE==0N) 

{ 

if  (ASCII==ON) 

{ 

for  (i=1;  i<=(BUFFER_TlME-cut)*F_SAMPLE;  i++) 
fprintf (fpHSE/'%e\n'',  HeanSqError[i3); 

>  /*  if  */ 

else  if  (BINARY==ON) 

{ 

if (fwrite( (char*) (MeanSqError+1), si zeof (float), (BUFFER  TIHE-cut)* 
F_SAMPLE,fpMSE)==(unsigned)(BUFFER_TIHE-cut)*F_SAMPLE)  ; 
else  if (ferror(fpMSE)  !=  0) 

ExitOnErrorC'Error  encountered  writing  error  data"); 
else  ExitOnErrorC'Unknown  error  handling  error  file."); 

>  /*  else  if  */ 

}  /*  if  */ 


if  (ASCII==0N) 

{ 

for  (i=1;i<=(BUFFER_TIME-cut)*F_SAMPLE;  i++) 
fprintf (fpOutSound,  "Xe\n",  outSOUNDCi]); 

>  /*  if  */ 

else  if  (B1NARY==0N) 

{ 

if (fwri ted  char*) (out SOUND+1), si zeof (float), (BUFFER  TlME-cut)* 
F_SAMPLE,fpOutSound)==(unsigned)(BUFFER_TIME-cut)*F_SAMPLE)  ; 
else  if (ferror(fpOutSound)  !=  0) 

ExitOnErrorC'Error  encountered  writing  output  acoustic  data"); 
else  ExitOnErrorC'Unknown  error  handling  acoustic  output  file."); 

>  /*  else  if  */ 
return(  0  ); 

/♦*★**  enO  putOutput  ***»*/ 

/»«t»*******>»*«r*«»««*«***  / 

/♦★*** 


*  FUNCTION;  dumpSpectrum( ) 

* 

*  This  function  handles  dumps  the 

* 

*  Arguments: 

* 

*  Return  value: 

* 

*  Functions  called: 

* 

*  Definitions  called: 

* 

* 

*  Global  variables  called: 

* 

*  Significant  memory  allocation: 

* 


signal  power  spectrum  (if  selected). 
pwrSpectrum 
0 

ExitOnErrorO 

ANSI  UNIX 

F_SAf1PLE  FFT_LENGTH 

none 

none 


Rif  defined  (  ANSI  ) 

int  dumpSpectrum  (  float  *pwrSpectrum  ) 
Relif  defined  (  UNIX  ) 
dumpSpectrum  (  pwrSpectrum  ) 
float  *pwrSpectrum; 

Rendif 

{ 

int  i; 

static  float  sequence=1 .0; 
char  f i leName[12],  mode[2]; 
static  FILE  *fp; 
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if Cf irstBuffer) 

{ 

if  (ASCII==ON) 
i 

■odeCO'W; 

•odeC13='f ; 

>  /*  if  */ 

else  if  (BINARY==ON) 

{ 

■odeC03= ' W ; 

■odeC1]='b'; 

>  /*  else  if  */ 

printf ("Enter  file  name  for  output  SPECTRUM  data: 

i f ( ( scanf ( "Xs" ,  f i leName) )==EOF ) 

ExitOnErrorC'Fatal  error  in  scanfO"); 
printf (“\n\n"); 

if  ((fp=fopen(fi leName,  mode))  ==  NULL) 

ExitOnError( "Error  opening  OUTPUT  SPECTRUM  data  file."); 

printf ("Select  output  format:  \n"); 
printfC  Enter  0  for  MATLAB  compatible  output. \n"); 
printfC  Enter  1  for  GRAFTOOL  compatible  output. \n"); 
if((scanf("Xi",  8i))==E0F) 

ExitOnErrorC'Fatal  error  in  scanfO"); 
if((i  1=  0)  S&  (i  !=  D)  ExitOnErrorCInvalid  output  selection") 
>  /*  if  */ 

if  <ASCII==0N) 

i f (f i rstBuf f er  S&  format==1) 

< 

fprintf (fp, "0  "); 

for  (i=1,  i<=FFT  LENGTH/2;  i++) 

fprintf(fp,"Xe  ",(f loat)(i*1)*F_SAMPLE/FFT_LENGTH); 
fprintf (fp,"\n"); 

}  /*  if  */ 

if  (format==1)  fprintf (fp,"%e  "/Sequence); 

for  (i=1;i<=FFT_LENGTH/2;  i++) 

fprintf (fp, "Xe  ",  pwrSpectrumCi]); 

fprintf (fp,"\n  "); 
sequence=sequence+1 .0; 

)  /♦  if  */ 

else  if  (BINART==ON) 

< 

if  (fwri  te(  (char*)  (puf-Spectrum+l  ),sizeof  (float  ),FFT_LENGTH/2,fp)= 
(unsigned) FFT_LENGTH/2)  ; 
else  if (ferror(fpOutSound)  !=  0) 

ExitOnErrorC'Error  encountered  writing  output  SPECTRUM  data") 
else  ExitOnError( "Unknown  error  handling  acoustic  SPECTRUM  file."); 

}  /♦  else  if  */ 
return(  0  ); 

/**♦*♦  end  dumpSpectrum  *♦♦♦*/ 

/liirtrkit 

*  FUNCTION:  processTiltO 

* 

*  This  function  handles  all  array  tilt  data.  It  calculates  the  X,  Y,  Z 

*  coordinates  of  each  hydrophone  as  a  function  of  time.  The  coordinate 

*  system  is  oriented  such  that  X  points  toward  the  signal  "origin"  and 

*  Z  points  down. 

* 
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*  Arguments: 

* 


*  Return  value: 

* 

*  Functions  called: 

* 

* 

*  Definitions  called: 

* 


ExitOnErrorO  MtrixO 
vectorO  free_MtrixO 

ANSI  UNIX 

DELTA_R  LOWER_SENSOR 

RADIAN  CHANNELS 

LOOK  DIRECTION  OFFSET 

tilt'ejffer 


*  Global  variables  called:  lastMinute 

* 

*  Significant  memory  allocation:  xx[][] 

*  2ZC]C] 

*  angleC] 

*  IdepthC] 


yy[3C] 

tiltC] 

udepthC] 


)»if  defined  (  ANSI  ) 

int  processTi Itlf loat  **x,  float  **y,  float  **2) 

#elif  defined  (  UNIX  ) 
processTi  lt(x,y,2'' 
float  **x,  **y ,  *Az; 

Oendif 

int  i,  j,  notEOF; 

float  *tilt,  *angle,  *udepth,  *ldepth,  **xx,  **yy,  theta 

char  f i  leNameC12]; 

FILE  *fp1,  *fp2; 

/*  Open  Data  Files  */ 

printf ("Enter  file  name  for  upper  tilt  sensor  data;  ")/ 

if  ( (scanf  ( "‘.s" ,  f  i  leName)  )==E0F ) 

ExitOnErrorC’Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ( <fp1=fopen(fi leName,  "rt"))  ==  NULL) 

ExitOnErrorC "Error  opening  UPPER  TILT  data  file."); 

if  (L0WER_SENS0R==0N) 

printf ("Enter  file  name  for  lower  tilt  sensor  data:  "); 

if ( (scanf ("%s" ,  f i leName) )==EOF) 

ExitOnErrorC  Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp2=fopen(fi leName,  "rt"))  ==  NULL) 

ExitOnErrorC'Error  opening  LOWER  TILT  data  file."); 
ldepth= ( float*) vector (TILT_BUFFER); 

>  /*  if  */ 

/*  Hemory  Allocation  */ 

t ilt= (float*) vector (TILT_BUFFER); 

ang I e= ( f loat* )vector(TI lT_BUF  F  ER ) ; 

udepth=(f loat*) vector (TILT_BUFFER); 

xx=(f loat**)matri x(CHANN£LS,  TILT_BUFFER); 

yy=(f loat**)matrix(CHANN£LS,  TILT_BUFFER); 

Z2=(f loat**)matrix(CHANNELS,  TILT_BUFFER); 


/*  read  upper  tilt  data  */ 


notEOF=1; 
whi le(notEOF) 

<. 

if (fscanf (fp1,"Xg  Xg  %g\n",8ti lt[i],6angle[i3,8udepthti])  !=  EOF)  i++; 
else  notE0F=0; 

> 


if  (LOWER  SENSOR==ON) 

{ 

j=1;  /*  read  lower  tilt  data  */ 

notE0F=1; 

whi le(notEOF) 

if(fscanf(fp2,  "%g\n’',8ldepth[j])  !=  EOF)  j++; 
else  notEOF=0; 

}  /*  whi le  */ 

if  {i<=j)  lastMinute=i-1; 

else  lastMinute=j-1; 

> 

else  lastHinute=i-1; 


/************ Th i s  is  the  assumed  array  geometry:  LINEAR***********/ 
for  (j=1;  j<=CHANNELS;  j++) 

for  (i=1;  i<=lastMinute;  i++) 

{ 


xx[j]Ci]=DELTA_R*(f  loat)(  3-1  )*cos(angle[:i  3/RADIAN)* 
sin(ti It [i 3 /RADIAN); 

yyCj3Ci3=DELTA_R*(f  loatX  j-1)*sin(angleCi3/RADIAN)* 
sinCti It[i3/RAD1AN); 

ZzCj3[:i3=DELTA  R*(float)(j-1)*cos(tiltCi3/RADIAN)+ 
OFFSET*cos(ti ltCi3/RADIAN)+udepth[i3; 

}  /*  for  */ 

>  /*  for  */ 


theta=(360.0-LOOK_DIR£CT10N)/RADIAN;  /*  coordinate  rotation  */ 
for  (3=1;  3<=CHANNELS;  3++) 

{ 

for  (i=1;  i<=lastMinute;  i++)  /*  points  x  into  signal  */ 

xC33Ci3=xxCi3[i3*cos(theta)-yy[33[i3*sin(theta); 

yCj3Ci3=xxCj3Ci3*sin(theta)+yyC33Ci3*cos<theta); 

zC33[i3=zz[33[i3; 

3  /*  for  */ 

}  /*  for  */ 


/*  Memory  Deallocation  */ 
if  (L0UER_SENS0R==0N)  freeCCchar*) Ldepth); 
free<(char*)ti It); 


> 

/■ 

* 

* 

* 

it 

* 


f reel (char* )angle); 
f  ree ((char*) udepth ) ; 
f ree_matrix(xx,  CHANNELS); 
f ree_matrix(yy,  CHANNELS); 
f ree_iiiatrix(zz,  CHANNELS); 
fclose(fpl); 

if  (LOUER_SENSOR==ON)  fclose(fp2); 
return(  0  ); 


/ 


END  process! i It 


FUNCTION:  processModesO 

This  function  handles  the  normal  mode  data.  It  calculates  hydrophone 
weights  and  group  speed.  The  user  must  insure  that  the  depth  vector 
and  eigenfunction  vector  are  of  equal  length. 
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* 


*  Arguments: 

zC][:] 

weight [] 

* 

ptrC 

*  Return  value: 

0 

*  Functions  called: 

vectorO 

ExitOnErrorO 

★ 

INTERPOLATE () 

dydxO 

*  Definitions  called: 

ANSI 

UNIX 

* 

PI 

SSP  LENGTH 

* 

EIGVAL_LENGTH 

ORDER 

* 

CHANNELS 

F_CARRIER 

*  Global  variables  called: 

Minute 

*  Significant  memory  allocation: 

depthC) 

Zntl 

* 

OnOff C] 

wCl 

* 

Kr[] 

dwdK 

* 

**★**/ 

#if  defined  (  ANSI  ) 

int  processModesC float  **z,  float  *weight,  float  *ptrC) 

#elif  defined  (  UNIX  ) 
processModesC z, weight, ptrC) 
float  **z,  *weight,  *ptrC; 

#endif 

< 

int  i,  j,  ptsEigVal,  set,  notEOF,  deadPhones,  weightNotAssigned; 
static  int  ptsEigFun; 
float  *w,  *Kr,  *dwdK,  err; 

static  float  depth[SSP_LENGTH+1],  Zn[SSP_LEN6TH+1],  OnOff [CHANNELS+1]; 
char  key,  f i leNameC12]; 

FILE  *fp1,  *fp2; 

if (firstBuffer==1) 

i 

w-(f loat*)vector(EI3VAL  LENGTH); 

Kr=(f loat*) vector (EIGVAL  LENGTH); 
dwdK= ( f loat*) vectorC  E IGVAL_LENGTH ) ; 

printf ("Enter  file  name  for  normal  mode  data  (eigenfunction):  "); 

if  ((scanf  ("/is",  f  i  leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp1=fopen(f i leName,  "rt"))  ==  NULL) 

ExitOnErrorC'Error  opening  NORMAL  MODE  data  file  (eigenfunction)") 

printf ("Enter  file  name  for  normal  mode  data  (eigenvalues);  "); 

if  ((scanf  ('"/.s",  f  i  LeName) )==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp2=fopen(fi leName,  "rt"))  ==  NULL) 

ExitOnErrorC'Error  opening  NORMAL  MODE  data  file  (eigenvalues).") 

i=1;  /*  read  normal  mode  data  */ 

notE0F=1; 

whi le(notEOF) 

if  (fscanf  (fpl ,  "/.g  %g  \n",  SdepthCi],  8Zn[i])  !=  EOF)  i++; 
else  notE0F=0; 

> 
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ptsEigFun=i-1; 

for(i=1; i<=ptsEigFun;i++)  depthCi]=depthCi3+CTD_0FFSET; 
i=1; 

notE0F=1; 
whi Le(notEOF) 

{ 

if(fsc8nf(fp2,  "%g  Xg  \n",  &wCi],  8KrCi3)  !=  EOF)  i++; 
else  notEOF=0; 

> 

ptsEigVal=i-1; 

for  (i=1;i<=CHANNELS;i++)  OnOff [i]=1 .0; 

printfC'Do  you  want  to  turn  off  any  hydrophones?  "); 

i f ( ( scanf ( "%s" , Skey ) )==E0F ) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

if(key=='y'  11  key=='Y') 

{ 

printf ("\nHow  many  hydrophones  must  be  secured?  "); 

i f ( (scanf ( "Xi " ,  &deadPhones) )==E0F ) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

for(i=1 ; i<=deadPhones; i++) 

printf {"\nEnter  hydrophone  number  to  secure:  ">; 

if((scanf("Xi",  &j))==E0F) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

if  (j>CHANNELS  11  j<1) 

ExitOnErrorC'Bad  hydrophone  identification"); 
OnOff Cj]=0.0; 

}  /*  for  */ 

)  /*  if  */ 

>  /*  if  */ 


for(i=1; i<=lastMinute; i++) 

{ 

if (depthCptsEigFun]<2r  CHANNELS][i]) 

{ 

printfC'Max  eigenfunction  depth  is:  Xf \n",depthCptsEigFun3); 
printfC'at  depthC]  index  of:  Xi\n",ptsEigFun); 
printfC'Max  depth  of  phone  number  Xi  is:  Xf\n\n", 

CHANNELS, zf  CHANNELS] [ i ] ) ; 

ExitOnErrorC'Fatal  data  set  error"); 

)  /*  if  */ 

>  /*  for  */ 

j=i; 

for  (i=1;  i<=CHANNELS;  i++) 

{ 

weightNotAssigned=1; 

while  (j<=ptsEigFun  88  weightNotAssigned) 

{ 

if (z[i]CMinute]<0.0  11  depth[j]<0.0) 

< 

printf("i=Xi\n",i); 

printf("j=Xi\n",p; 

printf ("Mi nute=Xi\n", Minute); 

printf ("zCi][Minute]  is:  XrVn",  zCiDCMinute]); 

printf ("depthCj]  is:  Xf \n",depth[ j]); 

printf ("Depth  less  than  zero  encountered  in  processModes . 

printf ("\n\n"); 
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ExitOnErrorC'Check  input  depths  for  coordinate  orientation") 
>  /*  if  */ 


if (depthCj]<z[i3CMinute]  &&  depth[j+1!]>zCi3CMinute]) 

{ 

set=>:0RDER+1)/2; 

lNTERP0LATE(8depthCj-set],  &ZnCj-set3,  OROER+1,  zCiDCMinute] 
&weight[i],  6err); 
wei ghtNot Ass igned=0; 

>  /*  if  */ 

)++; 

}  /*  whi le  */ 

>  /*  for  */ 

for  (i=1;  i<=CHANNELS;  i++)  weight[i3=0n0ff Ci3*weightCiD; 

if  (Minute==1) 

{ 

dydx (Kr , w, dwdK, pt sEi gVa  L ) ; 
for  (i=1; i<=ptsEigVal; i++) 

{ 

if  (uCi3<2.0*PI*F  CARRIER  &&  w[i+13>2.0*PI*f  CARRIER) 
set=(0RDER+1)/2; 

INTERPOLATE(Sw[i-set3,6dudK[i-set3,ORDER+1, 

2.0*PI*F  CARRIER, ptrC,&err); 

3  /*  if  */ 

3  /*  for  */ 


f ree( Cchar*)w); 
f  ree(  (ehar*)Kr); 
f ree((char*)dwdK); 

3  /*  if  */ 
returnC  0  ); 

/*n«r*  END  proeessdodes  ****♦/ 

/IhhhWt 

*  FUNCTION:  dydxO 

* 

*  This  function  estimates  derivatives. 

* 

* 
it 
it 
★ 

★ 

* 

* 

* 

* 

★ 

*  Global  variables  called: 

* 

*  Significant  memory  allocation: 


Arguments : 

Return  value: 
Functions  called: 
Definitions  called: 


xC3 

ddxC3 


ExitOnErrorO 

ANSI 

STEP 

none 

none 


y[3 

points 


UNIX 


#if  defined  (  ANSI  ) 

int  dydxCfloat  *x,  float  *y,  float  *ddx,  int  points) 

#elif  defined  (  UNIX  ) 

dydx(x,y,ddx, points) 

float  *x,  *y,  *ddx; 

int  points; 

#endif 

int  n; 
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tor  (n=1;n<=points;n++) 

{ 

if  ((n>=STEP)  &8  (n<=points-STEP))  /*center*/ 

ddx[n]={y[n+STEP]-yCn-STEP])/(xCn+STEP3-x[n-STEP3); 

else  if  (n<STEP)  /*beginning*/ 

ddx[n3=(yCn+STEP]-yC13)/(xCn+STEP]-xC13); 

else  if  (n>points-STEP)  /*end*/ 

ddxCn]=(y[points]-yCn-STEP])/(xCpoints]-xCn-STEPD); 

else 

ExitOnErrorC'lndex  error  in  dydx");  /*  sanity  check  */ 
>  /*  for  */ 
returnl  0  ); 

>  /SXXik*****'**********/ 

/*****  end  dydx  *****/ 

/***★* 

*  FUNCTION:  polintO 

A 

*  This  function  performs  polynomial  interpolation. 


★ 

*  Arguments: 

xaC] 

ya[3 

* 

n 

X 

* 

y 

dy 

"k 

*  Return  value: 

0 

★ 

*  Functions  called: 

vector () 

ExitOnErrorO 

* 

*  Definitions  called: 

ANSI 

UNIX 

* 

*  Global  variables  called: 

none 

* 

*  Significant  memory  allocation; 

dO 

cr3 

* 

k'kttkk/ 

Pif  defined  <  ANSI  ) 

int  polintC  float  *xa,  float  *ya. 

int  n,  float  x. 

float  *y,  float  *dy) 

Pelif  defined  (  UNIX  ) 
polint<xa,ya,n, x,y,dy) 
float  *xa,  *ya,  x,  *y,  *dy; 
int  n; 

Oendif 

{ 

int  i,  m,  ns=1; 

float  den,  dif,  dift,  ho,  hp,  w; 
float  *c,  *d; 

dif=fabs(x-xaC1 ] ); 

c=(f loat*) vector (n); 
d=(f loat*) vector  In); 

for  (i=1;  i<=n;  i++) 

{ 

if  ((dift=fabs(x-xaCi3))  <  dif) 

{ 

ns=i; 

dif=dift; 

>  /*  if  */ 

eCi]=yaCi]; 

dCi]=ya[i]; 

>  /*  for  ♦/ 

♦y^yaCns— ]; 

for  (m=1;  m<n;  m++) 
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for  (i=1;  i<=n-in;  i++) 
{ 


ho=x8[i]-x; 
hp=xaCi+B]-x; 
w=cCi+1]-d[i]; 
if  ((den=ho-hp)==0.0) 

ExitOnErrorC'Error  in  routine  POLINT"); 
den=w/den; 
d[i]=hp*den; 
cCi3=ho*den; 

>  /*  for  */ 


*y  +=  (*dy=(2*ns  <  (n-m)  ?  cCns+1]  :  dCns--])); 
>  /*  for  */ 
free((ehar*)d); 
free((char*)c); 
return (  0  ); 

/*****  end  polint  ****/ 

/***** 


*  FUNCTION;  window! ) 

* 

*  This  function  applies  a  Blackman 

* 

*  Arguments; 

* 

*  Return  value; 

* 

*  Functions  called; 

*  Definitions  called; 

* 

* 

*  Global  variables  called: 

* 

*  Significant  memory  allocation; 

* 


window  to  a  vector. 
data[]  N 

0 

none 

ANSI  UNIX 

PI 

none 

none 


*★***/ 


(Kif  defined  (  ANSI  ) 

int  window(float  *data,  int  N); 

#elif  defined  (  UNIX  ) 

window!  data,  N  ) 

float  *data; 

int  N; 

#endif 

{ 

int  n; 


for  !n=0;  n<N;  n++) 

dataCn+1]=dataCn+1]*!0.42+0.5*cos!2.0*PI*lf loat)!n-N/2)/!f loat)!N-1 )) 
+0.08*cos!4.0*PI*!float)!n-N/2)/lfloat)!N-1))>; 

>  /♦  for  */ 
return  !  0  ); 

/eeeee  end  w  i  ndow  eeeir*  / 


*  FUNCTION:  realft!) 

■k 

*  This  function  calculates  FFT's 

* 

*  Arguments:  dataC]  n 

*  isign 

* 

*  Return  value;  0 

* 
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* 

* 

Functions  called: 

fourl () 

* 

* 

Definitions  called: 

ANSI 

* 

* 

Global  variables  called: 

none 

* 

t 

Significant  aeaory  allocation: 

none 

***★*/ 

defined  (  ANSI  ) 

int  realft(float  *data,  int  n,  int  isign) 

delif  defined  (  UNIX  ) 

realftCdata,  n,  isign) 

float  *data; 

int  n,  isign; 

Aendif 

{ 

int  i,  i1,  i2,  i3,  i4,  n2p3; 

float  c1=0.5,  c2,  hlr.  Mi,  h2r,  h2i; 

double  wr,  wi,  wpr,  wpi,  uteap,  theta; 


> 

/■ 


theta=3.141592653589793/(double)n; 
if  (isign==1) 

{ 

c2  =  -0.5; 
fourKdata,  n,  1); 

)  /*  i<  */ 

else 

{ 

c2=0.5, 

theta  =  -theta; 

}  /*  else  */ 
wtemp-sin(0. 5*theta); 
wpr  =  -2.0*wtemp*wtemp; 
wpi=s'n(theta); 
wr=1 .  >wpr; 
wi*wp  , 

n2p3  2*n»3; 

for  '>1=2;  i<=ri/2;  i++) 

< 

i4=1  +  ( i3=n2p3-( i2=1  +  ( il^iri-l ))); 
Mr=c1*(data[i1]*dataCi3]); 
hi i=c1*(dataCi2]-d8ta[i4]); 
h2r  =  -c2'»(dataCi2]*data[i4]); 
h2i=c2*<dataCi1]-dataCi3]); 
dataC i1 ]=h1 r+wr*h2r-wi*h2i ; 
dataCi2]=h1 i*ur*h2i+wi*h2r; 
data[ i3]=h1 r-wr*h2r*wi*h2i; 
dataCiA]  =  -M  i+wr*h2irwi*h2r; 
wr=(wteinp=wr)*wpr-ui*wpi+«r; 
vi=wi*wpr+wte«ip*wpi+wi ; 

>  /  for  */ 
if  '  isign==1 ) 

dataC1]=(h1r=dataC1])+d8taC2]; 

dat8[2>h1r-dat8C2]; 

>  /♦  if  */ 


elsv 


detail ]= cl* ( (h1r=dataC1])roat8C2]); 

dataC2]=c1*(h1r-dataC2];; 

fourl (data,n,-1 ); 


}  /*  else  */ 

return  (  0  ); 


/**♦♦*  end  realft  *****! 


UNIX 
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*  FUNCTION:  fourK) 

* 

*  This  function  calculates  FFT's 


* 

*  Arguaents:  dataCl 

*  isign 

* 

*  Return  value;  0 

* 

*  Functions  called:  none 

* 

*  Definitions  called:  ANSI 

*  SWAP 

* 

*  Global  variables  called:  none 

* 

*  Significant  aemory  allocation:  none 

* 

Aif  defined  (  ANSI  ) 


int  fourKfloat  *data,  int  nn,  int  isign) 

Aelif  defined  (  UNIX  ) 

fourlldata,  nn,  isign) 

float  *data; 

int  nn,  isign; 

#endif 

int  n,  amax,  a,  j,  istep,  i; 

double  utemp,  wr,  wpr,  wpi ,  wi,  theta; 

float  tempr,  tempi. 


nn 


UNIX 


n=nn  «  1; 

for  (i=1;  i<n;  i+=2) 

< 

if  Ij  >  1) 

{ 

SUAPCdataC j],dataCi]); 
SUAP(dataC)+i],dataCi+1]); 

>  /*  for  */ 
m=n  »  1; 

whi  le  (m  >=  2  68  j  >  m) 

< 

j  -=  m; 

m  »=1; 

)  /*  whi le  */ 
j  I"; 

}  /*  for  */ 
mmax=2; 

whi le  (n  >  mmax) 

{ 

istsp=2*mmax; 

theta=6 . 2831 853071 7959 /< i s ign*mmax ) ; 

wtemp=sin(0. 5*theta); 

wpr  =  -2.0*wtemp*wtemp; 

wpi=sin(theta); 

wr=1 .0; 

wi=0.0; 

for  (ffi=1;  iiKmmax;  m+=2) 

{ 

for  (i=m;  i<=n;  i+=istep) 

{ 

j=i+mmax; 

te«pr=wr*data[]]-wi*dat8Cj+1]; 
tempi=wr*dataCj+1]+wi*data[j3; 
dataC j]=dataCi]-tempr; 
dataC j+1]=data[ 1+1]- tempi ; 
dataCi]  +=  tempr; 
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data[i+13  +=  tempi; 

}  /*  for  */ 

wr=(wtemp=wr)*wpr-ui*«pi+wr; 

wi=wi*wpr+wtemp*wpi+wi; 

>  /*  for  */ 
aMax=istep; 

)  /*  while  */ 
return  (  0  ); 

) 

/■ktirir*  end  foorl  *****/ 


*  FUNCTION;  vectorO 

* 

*  This  function  allocates  memory 

It 

*  Arguments: 

* 

*  Return  value; 
t 

it  Functions  called; 

* 

♦  Definitions  called: 

* 

*  Global  variables  called; 

* 

*  Significant  memory  allocation; 

* 


for  UNIT  OFFSET  vectors, 
length 

*v 

ExitOnErrorO 

ANSI  UNIX 

none 

vC] 


ftintirk/ 

P^1  defined  (  ANSI  ) 
float  *vector(int  length) 
Aelif  defined  (  UNIX  ) 
vectorC length) 
mt  length; 

Aendif 

{ 

float  *v; 


if  ((v=(f  loat»)mal  loci  ( length+1)*si2eof  (float  )))==NUU) 
ExitOnErrorC'Memory  allocation  failure  in  vectorO."); 
Aif  defined  (  ANSI  ) 


return  v; 

Aelif  defined  (  UNIX  ) 
return  (long  int)v; 


Aendif 


> 


/ 


END  vector  ****♦/ 
/*»* a****************** / 


*  FUNCTION:  matnxO 

♦ 

*  This  function  allocates  memory 

* 

*  Arguments; 

* 

*  Return  value; 

* 

*  Functions  called: 

A 

*  Definitions  called: 

* 

*  Global  variables  called: 

* 

*  Significant  memory  allocation: 

* 

Aif  defined  (  ANSI  ) 

float  AAmatrixlint  row,  mt  col) 


for  UNIT  OFFSET  a-D  arrays, 
row  col 

**m 

ExitOnErrorO 

ANSI  UNIX 

none 

m[]C] 
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#el.if  defined  (  UNIX  ) 

Htrix(rou,col.) 
int  row,  col; 

Hfendif 

int  i; 
float 

if  (((iF(f loat**)malloe( (unsigned) (row+1  )*si2eof (float*) ))==NULL) 
ExitOnErrorC'Allocation  failure  1  in  matrixO"); 
for  (i=1;  i<=row;  i++) 

{ 

if  ((mCi]=(f loat*)nia Hoc ( (unsigned) (col+1)*sizeof (float )))==NULL) 
ExitOnErrorC'Allocation  failure  2  in  aatrixO"); 

>  /*  for  */ 
iKif  defined  (  ANSI  ) 
return  m; 

IKelif  defined  (  UNIX  ) 
return  (long  int)m; 

#endif 

}  /»******«*******»**«««*/ 

/★****  end  matrix  *****/ 


*  FUNCTION:  f ree_raatrix() 

* 

This  function  deallocates  memory  from  UNIT  OFFSET  2-D  arrays. 


* 

* 

A 
* 

* 

* 

A 
A 
A 
A 
A 
A 
A 
A 

AAAAA/ 

#if  defined  (  ANSI  ) 

void  free_matrix(f loat  **m,  int  row) 

Aelif  defined  (  UNIX  ) 

f ree_matrix(m, row) 

float  *Am; 

int  row; 

Aendif 

{ 

int  i; 


Arguments:  mC]C] 

Return  value:  0 

Functions  called;  none 

Definitions  called:  ANSI 

Global  variables  called:  none 

Significant  memory  allocation:  none 


row 


UNIX 


> 


for(i=row;  i>=1;  i--) 
free((charA)m[i]); 
free((charA)m); 
return(  0  ); 


/AAAAA  end  free  matrix  aaaaa/ 


*  FUNCTION:  ExitOnError( ) 

A 

A  This  function  performs  an  abnormal  process  termination. 

A 

A  Arguments:  error_txtC] 

A 

A  Return  value:  none 

A 

*  Functions  called:  none 
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#elif  defined  (  UNIX  ) 
utrix(row,col) 

Int  row,  cot; 

Oendif 

int  i; 

float  **b; 


if  {<(iF(f loat**)BaUoc((unsigned)(row't-1)*si2eof (f loat*)))==NULL) 
ExitOnErrorC'Allocation  failure  1  in  aatrixO"); 
for  (i=1;  i<=row;  i++) 

{ 

if  ((m[i]=(f loat*)mal loc( (unsigned) (col+1)*sizeof(float)>)==NULL) 
ExitOnErrorC'Allocation  failure  2  in  aatrixO"); 

>  /*  for  */ 

•if  defined  (  ANSI  ) 


return  a; 

•elif  defined  (  UNIX  ) 
return  (long  intlm; 

•end if 

>  /**i>**j>*«iS**X*«ilA***«/>/ 

/*★***  end  matrix  ♦****/ 


/tkhkii 

*  FUNCTION:  f ree_matrix() 

* 

*  This  function  deallocates  memory  from  UNIT  OFFSET  2-D  arrays. 

* 

*  Arguments:  mC][]  row 


*  Return  value: 

* 

0 

*  Functions  called: 

* 

none 

t  Definitions  called; 

* 

ANSI 

*  Global  variables  called: 

* 

none 

*  Significant  memory  allocation 

* 

*****/ 

H'i  defined  (  ANSI  ) 

none 

void  free_matrix(f loat  **m,  int 
•elif  defined  (  UNIX  ) 
f ree_matrix(m, row) 
float  **m; 

row) 

int  row; 
•endif 


int  i; 


> 


for(i=row;  i>=1;  i — ) 
free((char*)mCi]); 
free((chap*)m); 
return(  0  ); 


/*★***  end  free  matrix  *****/ 
Z**************** «***«**»»** / 


*  FUNCTION:  ExitOnErrorO 

* 

*  This  function  performs  an  abnormal  process  termination. 

* 

*  Arguments:  error_txt[] 

* 

*  Return  value:  none 

* 

*  Functions  called:  none 
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it((c==COtiMA)  i;  (c==COLON))  c=SPACE; 
fputc(c,  out  fp); 

>  /*  while  */ 

}  /*  if  */ 

else  ExitOnErrorC’Error  opening  output  data  file"); 

>  /*  if  */ 

else  ExitOnErrorC'Error  opening  input  data  file"); 
fclose(sacB_fp); 
fclose<out_fp); 
exit(O); 

> 

l»if  (VERS10N==ANSI) 

void  ExitOnErrorlchar  error  txtC]) 

dlelif  (VERSION=»UNIX) 

ExitOnError(error_txt) 
char  error_txtC]; 

(It  end  if 
{ 

fprintfCstderr, "Program  run-time  error  ...\n‘'); 
fprintf (stderr,"Xs\n",error_txt); 
fprintf (stderr, " . . .now  exiting  to  system. .. \n"); 
exit(1); 

> 


C.  SACMl.C 

This  prograni  reads  the  data  output  by  the  upper  tilt 
sensor  and  formats  it  for  use  in  the  beamformer.  It 
calculates  60  second  averages  for  all  data  fields,  converts 
pressure  to  depth  and  calculates  tilt  direction  based  on 
current  velocity. 


/***★* 

*  PROGRAM  SACMl.C  vsn  1.0 

*  WRITTEN  BY;  Steven  Crocker 

*  LAST  UPDATE;  August  6,  1991 

* 

*  This  program  takes  SACM  data  from  the  Heard  Island  West  Coast 

*  Array  (upper  instrument  package)  and  condenses  it. 

* 

*  Pressure  is  converted  to  depth. 

* 

*  The  conductivity  is  not  processed.  All  values  output  are  60 

*  second  averages  of  the  input. 

* 


(fdef  ine 

VERSION  ANSI 

1* 

either  ANSI  or  UNIX 

*/ 

#def i ne 

SELECT  OUTPUT  OFF 

/* 

ON  or  OFF 

*/ 

#def ine 

C  25 

/* 

max  length  for  file  names 

*/ 

(Kdef  ine 

inBUFFER  1200 

1* 

input  buffer  size 

*/ 

c 

outBUFFER  50 

/* 

output  buffer  size 

*/ 

(fdefine  PI  3.14159265359 
(fdefine  RADIAN  57.2957795131 

(»include<stdio.h> 
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((include  <iMUoc.h> 

#include<Mth.h> 

(Kit  (VERSION—ANSI) 

(( inc  I  ude<std  I  i  b .  h> 

void  ExitOnErrorCehar  error_txtC]); 

#endif 

Mini) 

{ 

FILE  *saciii_fp,  *out_fp; 

char  c,  sacm_f i leCC3,  out_fileCC3; 

int  i,  j,  count,  loop_count,  inBUFFER_fuU,  out_ftag[6D; 

float  t1,  avgTilt,  avgPress,  avgVN,  avgVE,  avgTeap,  avgAngle, 
avgVelocity,  junk; 

static  float  inTi ItCinBUFFER],  inPressCinBUFFERl,  inVNCinBUFFER], 
inVECinBUFFER],  inTiaeCinBUFFER],  inTeapCinBUFFER], 
outTiaeCoutBUFFER],  outTi ItCoutBUFFERD, 
outPressEoutBUFFER],  outAngleCoutBUFFER], 
out Ve loci tyCoutBUFFER],  outTeopCoutBUFFER]; 

printf ("\n\n\nEntep  file  name  for  SACM  data.\n''); 
scant ("Xs",sacm_f i le); 
printf <"\n\n"); 

printf ("Enter  file  name  for  output.Xn"); 
scant ('7.s",out_fi  le); 

printf ("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n"); 
printfC'The  following  pertains  to  output  file  foraat.\n\n"); 

printf ("Units  of  measure:  Time  index  is  minutes. \n"); 
printfC  Angles  in  degreesXn"); 

printfC  Compass  direction  (TRUE)\n"); 

prinffC  Velocity  in  meters/second\n"); 

printfC  Temperature  in  degrees  C\n"); 

printfC  Depth  in  meters\n\n"); 

printf ("NOTE;  Output  column  order  is  same  as  above. \n\n\n"); 

for  (i=0;i<6;i++)  out_f lag[i]=1; 

#if  (SELECT_OUTPUT==ON) 

printf ("Include  TIME?  (y  or  n)\n"); 

scant ("%s", Sc); 

if(e==89  i;  c==121)  oot_f lagC0]=1; 
else  out_f l8gC0]=0; 
printf ("\n\n"); 

printf ("Include  TILT?  (y  or  n)\n"); 
scant  ("'/.s'^Sc); 

if(c==89  i;  c==121)  out_f lagC1]=1; 
else  out_f lagC1]=0; 
printf ("\n\n"); 

printf ("Include  CURRENT  DIRECTION?  (y  or  n)\n"); 
scant  CXs",  Sc); 

if(c==89  i;  c==121)  out_flagC2]=1; 
else  out_f lagC23=0; 
printf ("\n\n"); 

printfCInclude  CURRENT  VELOCITY?  (y  or  n)\n"); 
scant ("Xs", Sc); 

if(c==89  ;;  c==121)  out_flag[33=1; 
else  out_f lag[3]=0; 
printf("\n\n"); 

printfCInclude  TEMPERATURE?  (y  or  n)\n"); 
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scanf ("%s",8c); 

if(c==89  :i  c==121)  out_fLagC4]=1; 
else  out_f lagC4]=0; 
printf ("\n\n"); 

printf (“Include  DEPTH?  (y  or  n)\n"); 
scanf ("%s",6c); 

if(c==89  II  c==121)  out_flagC5]=1; 
else  out_f lag[5]=0; 

Kendif 

inBUFFER_full=inBUFFER; 

loop_count=0; 

if  ((sacm_fp=fopen(sacm_f 1 le,  "rt"))  !=  NULL)  ; 
else  ExitOnErrorC'Error  opening  input  data  file"); 
if  ((out_fp=fopen(out_file,  "wt"))  !=  NULL)  ; 
else  ExitOnErrorC'Error  opening  output  data  file"); 


whi  le 
< 


(1) 

for 

< 


(i=0;i<inBUFFER;i++) 


if  (fscanf (sacm_fp,  "%f  7,i  %f  7,f  Xi  %f  Zf\n“, 

8inTiiiie[i],&inVNCi],  SinVECi],  BinTempCiD, 
ginTiltri]^  fiinPressCi],  fijunk)  !=  EOF)  ; 
else  inBUFFER_full=i-2; 

>  /*  for  */ 
j=0; 

for  (i=0;i<inBUFFER  fuU;i++) 

{ 

t1=inTime[i]; 

avgTi lt=0.0; 

avgPress=0.0; 

avgVN=0.0; 

avgVE=0.0; 

avgTemp=0.0; 

eount=0; 

whi  le((inTiiiieCi]<t1+60.0)  &&  (i<inBUFFER  full)) 

{ 

avgTUt=avgTi  It  +  inTi  lt[i]/10.0; 
avgPress=avgPress  ♦  inPressCi]*1000.0; 
avgVN=avgVN  +  inVNC i ]/1000.0; 
avgVE=avgVE  +  inVE[i]/1000.0; 
av'gTemp=avgTeiiip  +  inTemp[i3/100.0; 
i++; 

count++; 

)  /*  whi le  */ 

avgT i lt=avgT i It /(f loat )count; 

avgPress=avgPress/(f  loat)count; 

avgTemp=avgT  emp/ (float) count ; 

avgVN=avgVN /(float) count ; 

avgVE=avgVE  /  ( f  loat )  count ; 

avgVe loc i ty=sqrt (avgVN*avgVN+avgVE*avgVE ) ; 

avgAngle=atan2(avgVE,  avgVN)*RAOIAN; 

if  (avgAngle<=0.0)  avgAngle=360.0+avgAngle; 

outTime[j]=(f loat) ( j+1 )+outBUFFER*loop_count; 

outTi lt[j)=avgTi It; 

outPressCj]=(avgPress-101 325. 0)/(9. 80665*1026.0); 

outTeiiip[i]=avgTemp; 

out Ve loc i tyCj ]=avgVe loci ty; 

outAngleC j]=avgAngle; 

)++; 

}  /♦  for  */ 


if  (inBUFFER_full  !=  inBurFER)  j  — ; 
for  ( 1=0; !<]; 1++) 

{ 

if  (out_f lagtO))  fpnntf (out_fp,  ">.6. Of",  outTiixeCi]) 
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> 


if  (out_f lagCI])  fprintf (out_fp, 
if  (out_f lagC2])  fpri-tf (out_fp, 
if  (out_f LagC3])  fprintf (out_fp, 
if  (out_f LagCA])  fprintf (out_fp, 
if  <out_f Lag[5])  fprintf (out_fp, 
fprintf (out_fp,  "\n"); 

}  /*  for  */ 


X8.5f",  outTiltCi]); 
%8.3f",  outAngleCi]); 
X8.6f"y  outVelocityCiD); 
%8.5f",  outTexpCiD); 
%8.0f",  outPress[i3); 


if  (inBUFFER_fulL  !=  inBUFFER)  exit(O); 
loop  count++; 

>  /*  while  */ 
exitd  ); 


#if  (VERSION==ANSn 

void  ExitOnErrorlchar  error_txt[]) 

Aelif  (VERSION==UNIX) 

ExitOnFrror(error_txt) 
char  error_txt[3; 

#endif 

fprintf (stderr/'Program  run-time  error  ...\n"); 
f print f( St derr, "%s\n", error_txt); 
fprintf  (stderr,". .  .now  exiting  to  system. .. Nn^'i^- 
exitCI); 

> 


D .  SACM2 . C 


This  program  takes  the  output  of  sacml.c  as  input.  It 
locates  and  copies  a  user  defined  subset  of  the  tilt  data  for 


use  in  beamforming. 


/♦**** 

*  PROGRAM  SACM2.C  vsn  1.0 

*  WRITTEN  BY;  Steven  Crocker 

*  LAST  UPDATE:  August  6,  1991 

* 

*  This  program  takes  data  processed  by  SACMI.C  and  cuts 

*  a  user  defined  segment  from  it.  The  segment  is  retained 

*  as  the  output  data  file.  The  input  file  is  not  affected. 

*  The  output  time  base  may  either  be  normalized  to  1  or  may 

*  retain  the  original  values. 

* 

*  The  option  is  given  to  accept  a  default  output  format.  This 

*  output  format  coincides  with  the  required  input  format  to  the 

*  time  domain  modal  beamformer. 

* 

*  The  output  elements  are  selectable  if  desired. 

♦****/ 

Adefine  VERSION  ANSI 
Adefine  C  25 
#include<stdio.h> 

^include  <malloc.h> 

«if  (VERSION==ANSI) 

Ainclude<stdl lb. h> 
void  ExitOnErrorCchar  error  txtC]); 

Pend  if 


/*  either  ANSI  or  UNIX  */ 

/*  max  length  for  file  names  */ 
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•ainO 

{ 

FILE  ♦input_fp,  *out_fp; 

char  c,  Default,  input_f  i  LeCC],  out_fUe[CD; 

int  i=0,  outi,  outN,  out_flag[73; 

float  Tilt,  Depth,  Time,  Temp,  Angle,  Velocity; 

printf ("\n\n\nEnter  file  name  for  input  data.\n”); 
scanf ("%s", input_f i le); 

if  ((input_fp=fopenCinput_file,  "rt"))  !=  NULL)  ; 
else  ExitOnErrorC'Error  opening  input  file"); 
printf ("\n\n"); 


printf ("Enter  file  name  for  output. \n"); 
scanf ("%s",out_fi le); 

if  ((out_fp=fopen(out_f i le,  "wt"))  !=  NULL)  ; 
else  ExitOnErrorC'Error  opening  output  file"); 


printf ("\n\n\nDo  you  want  the  default  output  conf iguration?\n"); 
printf ("This  option  formats  the  output  for  use  in  the  beamformer. \n") 
scanf ("Xs", Default); 

printf ("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n"); 
if(Default  !=  n'  &&  Default  !=  'N') 


printfC'The  following  pertains  to  output  file  format . \n\n"); 
printf ("Units  of  measure:  Time  Index  is  minutes. \n"); 
printfC  Angles  in  degreesXn"); 

printfC  Compass  direction  (TRUE)\n"); 

printfC  Velocity  in  meters/secondVn"); 

printfC  Temperature  in  degrees  C\n"); 

printfC  Depth  in  meters\n\n"); 

printf ("NOTE;  Output  column  order  is  same  as  above. \n\nVn"); 


printf ("Include  TIME?  (y  or  n)\n"); 
scanf  C/.s",c); 

if(c==e9  ;  c==121)  out_f  Lagt03=1; 
else  out_f lagC0]=0; 
printf  C\n\n"); 


printf ("Include  TILT?  (y  or  n)\n"); 
scanf  ("‘/.s",c); 

if(c==89  ;i  c==121)  out_flagC1]=1; 
else  out_f  lag[1]=0,■ 
prl^tf  C\n\n"); 

printfCInclude  CURRENT  DIRECTION?  (y  or  n)\n"); 
scanf  ( "',s" ,  c); 

if(c==89  ii  c==121)  out_f lagC23=1; 
else  out_f lagC2]=0; 
printf ("\n\n"); 

printfCInclude  CURRENT  SPEED?  (y  or  n)\n"); 
scanf ("Xs",c); 

if(c==89  II  c==121)  out_f lagC3]=1; 
else  out_f lag[3]=0; 
printf ("\n\n"); 


printfCInclude  TEMPERATURE?  <y  or  n)\n"); 
scanf  ("'/.s",c); 

if(c==89  II  c==121)  out_f lagC4]=1; 
else  out_f lag[4]=0; 
printf ("\n\n"); 


printfCInclude  DEPTH?  (y  or  n)\n"); 
scanf  (";.s",c); 

if(c==89  II  c^-=121)  out_f lag[5]=1; 
else  out_f lag[53=0, 

>  /*  if  */ 
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else 

{ 

out_f lag[0]=0; 
out_f  lagC13=1; 
out_f lag[23=1 ; 
out~ftagC33=0; 
out_f  lagC4D=0; 
out_f lag[53=1; 

>  /*  else  */ 

printf ("\n\nUhat  time  index  of  the  input  file  should  be  the  firstXn"); 
printf ("element  of  the  output  file?\n"); 
scanf ("%i",&out1); 
printf ("\n\n"); 

printf ("How  many  elements  should  the  output  file  contain?\n"); 
scanf ("X1",8outN); 
printf  ("\n\n"); 

if  (Default  !=  'n'  &&  Default  !=  'N') 

printf ("Should  the  time  base  be  normalized  to  begin  at  t=1?\n") 
scanf  ("/'.s",c); 

if(c==89  11  c==121)  out_f lagC63=1; 
else  out_f lag[63=0; 
printf ("\n\n\n"); 

3  /*  if  */ 

whi le  (1) 

( 

if  (fscanf (input_fp,  '7.f  %f  Xf  Xf  Xf  Xf  \n", 

8Time,  STilt/  &Angle,  ^Velocity,  STemp,  ftDepth)  !=  EOF) 

{ 

if((Time  >=  (float)outl)  &S  (Time  <  (f loatXoutl+outN))) 

{ 

if(out_flag[D]  &S  !out_flag[6J) 
fprintf (out_fp, "Xf ",  Time); 

else  if  (out_flagC03  &&  out_flag[63) 

fprintf (out_fp,  "Xf",  (f loat)(i+1)); 
if  (out_f lagC13)  fprintf (out_fp,  "  Xf",  Tilt); 

if  (out_f lagC23)  fprintf (out_fp,  "  Xf",  Angle); 

if  (out_f lagC33)  fprintf (out_fp,  "  Xf",  Velocity) 
if  (out_f lagC43)  fprintf (out_fp,  "  Xf",  Temp); 

if  (out_f lagC53)  fprintf (out_fp,  "  Xf",  Depth); 

fprintf (out_fp,  "\n"); 
i++  ; 

}  /*  if  */ 

else  if(Time  >=  (float ) (outl+outN) ) 

{ 

fclose(input_fp); 

fclose(out_fp); 

exit(O); 

>  /*  else  if  */ 

3  /*  if  */ 

else 

{ 

fclose(input_fp); 

fclose(out_fp); 

ExitOnErrorC'EOF  encountered  in  input  data  file"); 

3  /*  else  */ 

3  /*  whi le  */ 
exitd); 

3 

#if  (VERS10N==ANSI) 

void  ExitOnError(char  error  txtC3) 

Aelif  (VERSION==UNIX) 

ExitOnErrorCerror  txt) 
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char  error_txtC]; 

#endif 

{ 

fpr1ntf(stderr, "Program  run-time  error  ...\n"); 

fprintf(stderr/7.s\n",error_txt); 

fprintf (stderr,". . .now  exiting  to  system. .. \n"); 

exit(1); 

> 


E .  ARRATTEST . C 


This  program  performs  various  statistical  tests  on  each 
element  of  the  acoustic  array.  It  is  useful  for  isolating 
hydrophones  which  have  failed  during  the  experiment.  The 
statistical  tests  are  an  adaptation  of  those  found  in 
Reference  11. 


/tirkii 

*  PROGRAM:  ARRAYTST  usn  1.0 

*  WRITTEN  BY:  Steven  Crocker 

*  LAST  UPDATE:  October  21,  1991 

* 

It 


**♦*/ 


Adef ine 

UNIX  VERSION 

/* 

either  ANSI  or  UNIX 

*/ 

Adef ine 

CHANNELS  32 

/* 

number  of  hydrophones  on  array 

*/ 

Adef ine 

f  SAMPLE  228 

/* 

sampling  frequency 

*/ 

Adef ine 

BUFFER  TIME  60 

/* 

duration  of  time  averages 

*/ 

#include<stdio.h> 
#include<malloc.h> 
<!(incLude<math,  h> 


#if  defined  (  ANSI  ) 

#include<f loat .h> 

#include<stdl.ib.h> 
int  getlnput(void); 

int  momentsCf Loat  *data,  int  n,  float  *ave,  float  *ave2,  float  *adev, 
float  *sdev,  float  *svar,  float  *skew,  float  *curt) 
float  *vector(int  length); 
float  **matrix(int  row,  int  col); 
int  f ree_matrix(f loat  **m,  int  row); 
void  ExitOnErrorlchar  error_txt[]); 

Aelif  defined  (  UNIX  ) 
int  vectorO, 

Mtnxl ), 
get  Input  O, 

■omentsO, 
f ree_matrix(), 

ExitOnErrorC ); 

Aendif 

/♦  Global  Variable  Declarations  */ 
float  AAinSOUND; 
int  TirstBuffer=1; 

mainO 

{ 

int  1,  j=1,  k,  flagCSl; 
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float  ave,  ave2,  adev,  sdev,  svar,  skew,  curt,  *te«p; 
char  f i leNaneCSO],  key; 

FILE  *fp1,  *fp2,  *fp3,  *fp4,  *fp5,  *fp6,  *fp7; 

inSOUND=(f loat  **)aatrix(CHANNELS,  BUFFER  TIHE*F  SAMPLE); 
teap=(f loat  *)vector(BUFFER_TIME*F_SAMPLET; 

for  (i=0;i<=7;i++)  flag[1>0; 

printf ("Answer  y  or  n  to  the  following  output  options. \n"); 
printf ("Average  value  (  over  %i  seconds  );  ",  BUFFER_TIME); 
scanfC'Zs",  &key); 
if(key==89  |;  key==121) 

{ 

flag[1>1; 

printf ("Enter  file  name;  "); 
scanfC'Zs",  fileName); 
printf ("\n\n"); 

if  ((fp1=fopen(f i leName,  "wt"))  ==  NULL) 

ExitOnError( "Error  opening  file."); 

>  /*  if  */ 

printfC'Mean  squared  value  (  over  Zi  seconds  ):  ",  BUFFER_T1ME); 
scanfC'Zs",  8key); 
if(key==89  J)  key==121) 

{ 

f lagC2]=1; 

printf ("Enter  file  name:  "); 
scanfC'Zs",  fileName), 
printf("\n\n"); 

if  ((fp2=fopen(f i leName,  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  file."); 

}  /*  if  */ 

printf ("Average  deviation  (  over  Zi  seconds  ):  ",  BUFFER_TIME); 
scanfC'Zs",  &key); 
if(key==89  key==121) 

flagC3]=1; 

printf ("Enter  file  name:  "); 
scanfC'Zs",  fileName); 
printf ("\n\n"); 

if  ((fp3=fopen(fi leName,  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  file."); 

)  /*  if  */ 

printf ("Standard  deviation  (  over  Zi  seconds  );  ",  BUFFER_TIME); 
scanfC'Zs",  Rkey); 
if(key==89  H  key==121) 

{ 

flagC4]=1; 

printf ( "Enter  file  name:  "); 
scanfC'Zs",  fileName); 
printf ("\n\n"); 

if  ((fp4=fopen(fi leName,  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  file."); 

}  /*  if  */ 

printf ("Variance  (  over  Zi  seconds  ):  ",  BUFFER_TIHE); 
scanfC'Zs",  &key); 
if(key==89  li  key==121) 

{ 

fl8g[5]=1; 

printf ("Enter  file  name:  "); 
scanfC'Zs",  fileName); 
printf <"\n\n"); 

if  ((fp5=fopen(fi leName,  "wt"))  ==  NULL) 

ExitOnErrorC’Error  opening  file.”); 

>  /*  if  */ 
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printf ("Skewness  (  over  Xi  seconds  ):  ",  BUFFER_T1HE); 
scant ("Xs",  ikey); 
if(key==89  1|  key==121) 

{ 

flag[6]=1; 

printf ("Enter  file  nane: 

scanfC'Xs",  fileNane); 

printf ("\n\n"); 

if  ((fp6=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 


printf ("Kurtosis  (  over  Xi  seconds  ):  ",  BUFFER_T1ME); 
scanfC'Xs",  &key); 
if(key==89  ||  key==121) 

flagC7]=1; 

printf ("Enter  file  name:  "); 
scanfC'Xs",  fileName); 
printf ("\n\n"); 

if  ( (f p7=fopen(f i leName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 


while  (j<120) 

{ 

getInputO; 

if  (flagCI])  fprintf(fp1,"Xi", j); 
if  (flagC2])  fprintf (fp2, "Xi", j); 
if  (flagC3])  fprintf (fp3,"Xi",j); 
if  (flagC4])  fprintf(fp4,"Xi", j); 
if  (flagCSD)  fprintf(fp5,"Xi", j); 
if  (flagC6])  fprintf(fp6,"Xi", j); 
if  (flagC7])  fprintf (fp7,"Xi",j); 


for  (i=1;i<=CHANNELS;i++) 

< 

for(k=1;k<=BUFFER_nME*F_SAMPLE;k++) 
tempt k]=inSOUNDCi]Ck]; 


moment s( temp, BUFFER_TIME*F_SAMPLE,&ave,&ave2, 
6adev,Ssdev,&svar,&skew,&curt); 


if 

(flagtl])  fprintf(fp1. 

"  Xf 

",  ave); 

if 

(flag[2])  fprintf(fp2. 

"  Xf 

",  ave2); 

if 

(flagC3])  fprintf(fp3. 

"  Xf 

",  adev); 

if 

(flagC4])  fprintf(fp4. 

"  Xf 

",  sdev); 

if 

(flagCS])  fprintf(fp5, 

"  Xf 

",  svar); 

if 

(flag[6])  fprintf(fp6. 

"  Xf 

",  skew); 

if 

(flagC7])  fprintf(fp7. 

"  Xf 

",  curt); 

)  /*  for  */ 


if  (■'.iqtl])  fprintf (fp1,"\n"); 
if  ■fi,f.’gC2])  fprintf(fp2,"\n"); 
if  (fi.agC33)  fprintf {fp3,"\n''); 
if  (flagCA])  fprintf (fp4,"\n"); 
if  (flagCS])  fprintf (fp5,"\n"); 
if  (flagC6))  fprintf (fp6, "\n"); 
if  (flagC7])  fprintf (fp7, "\n"); 
f irstBuffer=0; 

)**; 

>  /*  whi le  */ 


exit(O); 

END  main  ****«»*»********AAA>AtAAAAA«>t/ 


*  FUNCTION:  get  Input () 
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* 

*  This  function  handles  all  acoustic  input.  It  also  provides  one  of 

*  two  nornal  process  terminations  available  in  the  program. 


* 

*  Arguments: 

it 

none 

*  Return  value; 

0 

A 

*  Functions  called: 

vector!) 

ExitOnErrorO 

m 

t  Definitions  called: 

ANSI 

UNIX 

t 

F  SAMPLE 

CHANNELS 

it 

BUFFER_TIME 

t 

*  Global  variables  called; 

inSOUNDCDC] 

f irstBuffer 

A 

*  Significant  memory  allocation: 

diskBufferCi 

* 

Uii  defined  (  ANSI  ) 
int  getlnput(void) 

#elif  defined  (  UNIX  ) 
get  Input () 

Aendif 

< 

int  i,  ],  items,  buffer; 
float  AdiskBuffer; 
char  f i leNameCSO]; 
static  f I'  E  AfpStatic; 

if  (f irstBuffer) 

{ 

pnntf ("Enter  file  name  for  input  acoustic  data;  "); 

if ( (scanf ( "^s",  f i leName) )==E0F ) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ( (fpStat ie=fopen(fi leName,  "rb"))  ==  NULL) 

ExitOnErrorCError  opening  INPUT  ACOUSTIC  data  file."); 

>  /*  if  */ 

/*  Memory  Allocation  */ 

buffer=BUFFER_TIME*F_SAMPLE*CHANNELS; 
diskBuffer=(f loat*) vector (buffer); 

1 tems=f read( ( char*) (di skBuf f er+1 ) ,  siteof (float),  buffer,  fpstatic); 
if  ( items==buffe'')  /*continue*/; 
else  if (ferror(fpStatic)  !=  0) 

Ex  1 tOnError ( "Error  encountered  while  reading  input  acoustic  date") 
else  if (feof (fpStatic)  !=  0) 

printf("\n\t**««*  **««**'***************************««***** '««***\n") 

printfC'Vt  End  of  File  reached:  EXECUTION  COMPLETEXn"); 
pnntf  ("Nt**************************************************)^''^) 

fclose(fpStatic); 

exit(O); 

}  /*  else  if  */ 

else  ExitOnEnrorC'Unknown  error  handling  acoustic  input  file."); 

for  (1=1;  i<=BUFFER  TIHE*F  SAMPLE;  i++) 

{ 

for  (j=1;  j<=CHANNELS;  j++) 

inSOUNOCjlCi]  =  disk3ufferCCHANNELS*(i-1)+3]; 

>  /*  for  */ 

/*  Deallocate  Memory  */ 
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free(<char*)d1skBuffer); 
return(  0  ); 

>  XU******** *«****»»***»*«*/ 

/*****  end  get  Input  ***★*/ 
/to*********************/ 


*  FUNCTION:  aoiientsO 

* 

*  This  function  calculates  the  aean,  aean  squared  value,  average 


*  deviation,  standard  deviati 

*  of  a  data  vector. 

* 

*  Arguments; 

* 

* 

* 

* 

* 

*  Return  value: 

* 

*  Functions  called: 

* 

*  Definitions  called; 

A 

*  Global  variables  called: 

* 

*  Significant  memory  allocation 

* 


,  variance,  skewness  and  kurtosis 


data  n 

ave  ave2 

adev  sdev 

svar  skew 

curt 

0 

ExitOnErrorO 

ANSI  UNIX 

none 

none 


#if  defined  (  ANSI  ) 

int  momentslf loat  *data,  int  n,  float  *ave,  float  *ave2,  float  *adev, 
float  *sdev,  float  *svar,  float  *skew,  float  *curt) 

Aelif  defined  <  UNIX  ) 

moment s (data, n, ave, ave2, adev, sdev, svar, skew, curt) 
int  n; 

float  *data,  *ave,  *ave2,  *adev,  *sdev,  *svar,  *skew,  *curt; 

Aend i f 
{ 

int  j; 
float  s,p; 

if  (n<=1)  ExitOnErrorC'n  must  be  at  least  2  in  momentO"); 
s=0.0; 

*8ve2=0.0; 

for  ( i='l;  i<=n;  j++) 

{ 

s  +=  dataCj]; 

*ave2  +=  data[j]*data[j]; 

}  /*  for  */ 

*8ve=s/n; 

*ave2  /=  n; 

*adev=(*svar)=(*skew)=(Acurt )=0.0; 
for  <j=1; j<=n,- j++) 

{ 

♦adev  +=  fabs(s=dataC j ]-(*ave) ); 

♦svar  +=  (p=s^s); 

♦skew  +=  <p  ♦=  s); 

♦curt  +=  (p  ♦=  s); 

>  /♦  for  ♦/ 

♦adev  /=  n; 

♦svar  /=  (n-1); 

♦sdev=sqrt (♦svar ) ; 
if  (♦svar) 

{ 

♦skew  /=  (ne(^svar)^(^sdev)); 

♦curt= (♦curt )/(n^ (♦svar )^(^svar) ) -3.0; 

}  /♦  if  ♦/ 

else  ExitOnErrorC'No  skew/kurtosis  when  variance  =  0  (in  momentO)") 
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return(  0  ); 

>  />*««A**»»***«**A******/ 

/**★**  end  aonent  *****/ 
/aaaaaaaaaaaaaaaaaaaaa*/ 

/**★** 

A  FUNCTION;  vector () 

* 

*  This  function  allocates  aeaory 

* 

*  Arguaents: 

* 

*  Return  value: 

* 

*  Functions  called: 

* 

*  Definitions  called: 

* 

*  Global  variables  called; 

* 

*  Significant  aeaory  allocation; 

* 


for  UNIT  OFFSET  vectors, 
length 

*v 

ExitOnErrorO 

ANSI  UNIX 

none 

vC] 


#if  defined  (  ANSI  ) 
float  *veetor(int  length) 
Aelif  defined  (  UNIX  ) 
vector! length) 
int  length; 

Aendif 

< 

float  ; 


if  ((v=(f  loat*)maUoc((  Length+1  )*sizeof  (float  )))==NULL) 
ExitOnErrorC "Memory  allocation  failure  in  vectorO."); 
Aif  defined  (  ANSI  ) 
return  v; 

Aelif  defined  (  UNIX  ) 
return  (long  int)v; 

Aendif 

}  /AAAA-*A*A-AAAAAAAAAAAAA*/ 

/tirlrkii  END  vector  **★♦*/ 


/***♦* 

*  FUNCTION:  matrix!) 

* 

*  This  function  allocates  memory 

* 

*  Arguments; 

A 

*  Return  value: 

* 

*  Functions  called: 

* 

*  Definitions  called: 

* 

*  Global  variables  called: 

* 

*  Significant  memory  allocation: 

* 


for  UNIT  OFFSET  2-D  arrays, 
row  col 

e*m 

ExitOnErrorO 

ANSI  UNIX 

none 

mC]C] 


Aif  defined  (  ANSI  ) 

float  **matrix(int  row,  int  col) 

Aelif  defined  (  UNIX  ) 

matrix(row,col) 

int  row,  col; 

Aendif 

int  i; 

float  **m; 
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if  ((iFCf loat**)ma I loc( (unsigned) (row+1)*sizeof (float*) ))==NULL) 
ExitOnErrorC'AUocation  failure  1  in  aatrixO"); 
for  (i=1;  i<=row;  i++) 

{ 

if  ((BCi]=(f loat*)i»a I loc( (unsigned) (col+1)*sizeof (float )))==NULL) 
ExitOnErrorC'AUocation  failure  2  in  BatrixO"); 

>  /*  for  */ 

Ifif  defined  (  ANSI  ) 
return  m; 

#elif  defined  (  UNIX  ) 
return  (long  int)m; 

#endif 

>  />****« AAA* AAA* A***/ 

/*****  end  matrix  *★***/ 

/AAAAAAAAAAAAAAAAAAAAAA/ 


*  FUNCTION:  free  r»atrix() 

A 

*  This  function  deallocates  memory  from  UNIT  OFFSET  2-D  arrays. 


A 

*  Arguments: 

mC]C] 

rc 

A 

*  Return  value: 

0 

A 

*  Functions  called: 

none 

A 

A  Definitions  called: 

ANSI 

UNIX 

A 

*  Global  variables  called; 

none 

A 

*  Significant  memory  allocation 

none 

A 

AAA**/ 

#if  defined  (  ANSI  ) 

void  free_matrix(f loat  **m,  int 

row) 

#elif  defined  (  UNIX  ) 

f  ree_matrix(in,  row) 
float  **m; 
int  row; 

Aendif 

< 

int  i; 


> 


for(i=row;  i>=1;  i — ) 
f reel (char*)m[ i ]  ); 
f reel (ch8r*)m); 
returnl  0  ); 


/ 


/****•  END  free  matrix  ***** / 


*  FUNCTION:  ExitOnErrorO 

* 

*  This  function  performs  an  abnormal  process  termination. 


*  Arguments:  error  txtC] 

A 

*  Return  value:  none 

A 

*  Functions  called:  none 

A 

*  Definitions  called;  ANSI 

A 

A  Global  variables  called:  none 

A 

A  Significant  memory  allocation:  none 

A 

/ 


UNIX 


85 


#if  defined  (  ANSI  ) 

void  ExitOnError(char  error_txtC]) 

«elif  defined  (  UNIX  ) 

ExitOnError(error_txt) 
char  error_txtC]; 

#endif  * 

{ 

fprintflstderr, "Program  run-time  error  ...\n"); 
f pri ntf ( stderr, "%s\n",error_txt); 

fprintf (stderr,". . .now  exiting  to  system. . .\n");  , 

exit(O); 

}  /eexossxxixxxxxxexxxxxA*****/ 

/*★***  END  ExitOnError  *****/ 

/AAAAAAAAAAAAAAAX***********/ 


) 
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